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Abstract. This review demonstrates that neutrino emission from dense hadronic 
component in neutron stars is subject of strong modifications due to collective 
effects in the nuclear matter. With the most important in-medium processes incor- 
porated in the cooling code an overall agreement with available soft X ray data can 
(f) ■ be easily achieved. With these findings so called "standard" and "non-standard" 

CN| \ cooling scenarios are replaced by one general "nuclear medium cooling scenario" 

which relates slow and rapid neutron star coolings to the star masses (interior den- 
C*) ' sities). In-medium effects take important part also at early hot stage of neutron 

star evolution decreasing the neutrino opacity for less massive and increasing for 
\m~) ' more massive neutron stars. A formalism for calculation of neutrino radiation from 

^sL nuclear matter is presented that treats on equal footing one-nucleon and multiple- 

nucleon processes as well as reactions with resonance bosons and condensates. 

o 
o 

O ■ 1 Introduction 

O ,■ The EINSTEIN, EXOSAT and ROSAT observatories have measured surface 

temperatures of certain neutron stars (NS) and put upper limits on the sur- 
face temperatures of some other NS ( cf. [1-3] and further references therein). 

£i | The data for some supernova remnants indicate rather slow cooling, while the 

d • data for several pulsars point to an essentially more rapid cooling. 

Physics of NS cooling is based on a number of ingredients, among which 
• ■ *h . the neutrino emissivity of the high density hadronic matter in the star core 

K% ' plays a crucial role. Neutron star temperatures are such that, except first 

i3 . minutes-hours, neutrinos/antineutrinos radiate energy directly from the star 

without subsequent collisions, since A„,A P 3> -R, where A„,Ap are the neu- 
trino and antineutrino mean free paths and R is star radius. In the so called 
"standard scenario" of the NS cooling (scenario for slow cooling) the most im- 
portant channel up to temperatures T ~ 10 9 K belongs to the modified Urea 
(MU) process nn — > npeD. First estimates of its emissivity were done in 
[4,5]. References [6,7] recalculated the emissivity of this process in the model, 
where the nucleon-nucleon (NN) interaction was treated with the help of 
slightly modified free one-pion exchange (FOPE). This important result for 
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the emissivity, e v [FOPE] , was proved to be by an order of magnitude larger 
than the previously obtained one. Namely the value e, y [FOPE] was used in 
various computer simulations resulting in the "standard scenario" of the cool- 
ing, e.g. cf. [8-10]. Besides the MU process, in the framework of the "standard 
scenario " numerical codes included also processes of the nucleon (neutron and 
proton) bremsstrahlung (NB) nn — > nnvv and np — > npvD, which lead to a 
smaller contribution to the emissivity then the MU, cf. [11,6]. Medium effects 
enter the above two-nucleon (MU and NB) rates mainly through the effec- 
tive mass of the nucleons which has a smooth density dependence. Therefore 
within FOPE model the density dependence of the reaction rates is rather 
weak and the neutrino radiation from a NS depends only very weakly on its 
mass. This is the reason why the "standard scenario" based on the result [6], 
though complying well with several slowly cooling pulsars, fails to explain the 
data of the more rapidly cooling ones. Also "standard scenario " included pro- 
cesses contributing to the emissivity in the NS crust which become important 
at a lower temperature. 

The non-standard scenario included so called exotica, associated with 
different types of direct Urea-like processes, i.e. the pion Urea (PU) [12] 
and kaon Urea (KU) [13,14] /3-decay processes and direct Urea (DU) on 
nucleons and hyperons [15] possible only in sufficiently dense interiors of 
rather massive NS. The main difference in the cooling efficiency driven by the 
DU-like processes on one hand and the MU and NB processes on the other 
hand lies in the rather different phase spaces associated with these reactions. 
In the MU and NB case the available phase space is that of a two-fermion 
origin, while in the pion (kaon) /3-decay and DU on nucleons and hyperons 
it is that of a one-fermion origin. Critical density of pion condensation in NS 
matter is g C7r ~ (1 -j- 3) go depending on the type of condensation (neutral or 
charged) and the model, sec [16-19]. Critical density of kaon condensation 
is q c k — (2 -7- 6)go depending on the type (K~ or K°, S or P wave) and 
the model, see [20,21]. Critical density for the DU process is g c u ~ (2 -=- 6)go 
depending on the model for the equation of state (EQS), see [15,18]. Recent 
calculations [19] estimated critical density of neutral pion condensation as 
2.5go and for the charged one as 1.7 go, whereas variational calculations [18] 
argued even for smaller critical densities (~ 1.3 go for 7r° condensation). On 
the other hand the EQS of [18] allows for DU process only at g > 5go- 

There is no bridge between "standard" and "non-standard" scenarios due 
to complete ignorance of in-medium modifications of TV TV interaction which 
allows for strong polarization of the soft modes (like virtual dressed pion 
and kaon modes serving a part of in- medium baryon-baryon interaction). 
Only due to enhancement of such a polarization with the baryon density 
mentioned condensates may appear and it seems thereby quite inconsistent 
to ignore these softening effects for g < g C7Tl g c x, and suddenly switch on the 
condensates for g > g C7Tl g c x- 



Medium Effects 3 

Now let us basing on the results [22-25,16,26] briefly discuss a general 
"nuclear medium cooling scenario " which treats obvious caveats of two men- 
tioned above scenarios. First of all one observes [22] that in the nuclear mat- 
ter many new reaction channels arc opened up compared to the vacuum 
processes. Standard Feynman technique fails to calculate in-medium reac- 
tion rates if the particle widths are important since there are no free particle 
asymptotic states in matter. Then summation of all perturbative Feynman 
diagrams where free Green functions are replaced by the in-medium ones 
leads to a double counting due to multiple repetitions of some processes (for 
an extensive discussion of this defect see [27]). This calls a formalism dealing 
with closed diagrams (integrated over all possible in-medium particle states) 
with full non-equilibrium Green functions. Such a formalism was elaborated 
in [23,24] first within quasiparticle approximation (QPA) for nucleons and was 
called in [24] "optic theorem formalism (OTFj in non- equilibrium diagram 
technique". It was demonstrated that standard calculation of the rates via 
squared reaction matrix elements and calculation using OTF coincide within 
QPA picture for the fermions. In [28] the formalism was generalized to in- 
clude arbitrary particle widths effects. The latter formalism treats on equal 
footing onc-nucleon and multiplc-nucleon processes as well as resonance re- 
action contributions of the boson origin, as processes with participation of 
zero sounds and reactions on the boson condensates. Each diagram in the 
series with full Green functions is free from the infrared divergences. Both, 
the correct quasiparticle (QP) and quasiclassical limits are recovered. 

Except for very early stage of NS evolution (minutes - hours) typical 
averaged lepton energy (> T) is larger then the nuclcon particle width /V ~ 
T 2 /efn and the nucleons can be treated within the QPA. This observation 
much simplifies consideration since one can use an intuitive way of separation 
of the processes according to the available phase space. The one-nucleon 
processes have the largest emissivity (if they are not forbidden by energy- 
momentum conservations), then two-nucleon processes come into play, etc. 

In the temperature interval T c < T < T opac (T c is typical tempera- 
ture for the nucleon pairing and T opac is typical temperature at which neu- 
trino/antineutrino mean free path A^/Ap is approximately equal to the star 
radius R) the neutrino emission is dominated by the medium modified Urea 
(MMU) and medium nucleon bremsstrahlung (MNB) processes if one-nucleon 
reactions like DU, PU and KU are forbidden, as it is the case for g < 
QcUi Qcix, QcK- Corresponding diagrams for MMU process are schematically 
shown in Fig. 1. References [22-25,16] considered NN interaction within 
Fermi liquid Landau-Migdal approach. They incorporated the softening of 
the medium one-pion exchange (MOPE) mode, other medium polarization 
effects, like nucleon-nucleon correlations in the vertices, renormalization of 
the local part of NN interaction by the loops, as well as the possibility of 
the neutrino emission from the intermediate reaction states and resonance 
DU-like reactions going on zero sounds and the boson condensates. Refer- 
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Fig. 1. Antineutrino emission from a nucleon leg (left graph) and from intermediate 
scattering states (right) in MMU process. Full dot includes weak coupling vertex 
renormalization. 



ences [23,25,16] have demonstrated that for g > qq second diagram of Fig. 
1 gives the main contribution to the emissivity of MMU process rather than 
the first one which contribution has been earlier evaluated in the framework 
of FOPE model in [6] . This fact essentially modifies the absolute value as well 
as the density dependence of the nn — > npev process rate which becomes to 
be very sharp. Thereby, for stars of masses larger than the solar mass the 
resulting cmissivities were proved to be substantially larger than those values 
calculated in FOPE model. With increase of the star mass (central density) 
pion mode continues to soften and MMU and MNB rates still increase. At 
Q > Qc-k pion condensation begins to contribute. Actually, the condensate 
droplets exist already at a smaller density in the mixed phase appearing in 
systems having more than one conserved charge [29,30]. At T > T me it, where 
Tmeit is the melting temperature, roughly ~ several MeV, the mixed phase 
is in liquid state and PU processes on independent condensate droplets are 
possible. At T < T me i t condensate droplets are placed in a crystalline lattice 
that substantially suppresses corresponding neutrino processes. 

Reference [12] considered the reaction channel n — ► pTr~ev, whereas [22,23] 
included other possible pion ir + , n , ir° condensate processes with charged 
and neutral currents (e.g., like nir® — > pev, nir^ — » pvv and nir® — > nvv) 
as well as resonance reactions going on zero sounds which are also possible 
at q < g c1T . Due to NN correlations all pion condensate rates are signifi- 
cantly suppressed (by factors ~ 10 - 100 compared to first estimate [12], see 
[31,22,23,32]. At g ~ g c1T both MMU and PU processes are of the same order 
of magnitude [23] demonstrating a smooth transition to higher densities (star 
masses) being absent in the "standard" and "non-standard" scenarios. 

For T < T c the reactions of neutrino pair radiation from supcrfhiid 
nucleon pair breaking and formation (NPBF) shown in Fig. 2 become to 
be dominant processes. The neutron pair breaking and formation (nPBF) 
process for the case of the lSo pairing was first calculated in [33] using 
standard Bogolyubov technique. Later this process was independently cal- 
culated in [24,25] as demonstration of efficiency of OTF within the closed 
non-equilibrium diagram technique developed there. Moreover [24] calculated 
emissivity of the corresponding process on proton (pPFB) taking into account 
strong coupling pvv vertex renormalization (see first diagram (17) below). 
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Fig. 2. Neutrino-antineutrino emission from Cooper pair-breaking (left graph) and 
pair- formation processes (right graph). 



It results in one-two order of magnitude enhancement of pPFB cmissivity 
compared to that would be estimated with the vacuum vertex, leading to 
that both nPFB and pPFB emissivities can be of the same order of mag- 
nitude. Emissivities of NPBF processes have the same suppression factor 
~ exp(-2 A/T) as MIX, NB, MMU and MNB at T < T c but compared to 
the latter the NPBF processes have a large one nucleon phase space volume. 
References [34,35] included peculiarities of 3P2 pairing. In the 3P2G m J | = 2) 
case, where mj is projection of the total pair momentum onto quantization 
axis, exponential suppression of the specific heat and the emissivity is re- 
placed by only a power law suppression since the gap vanishes at the Fermi 
sphere poles, the possibility first remarked in [24] and then in [16]. 

At this instance I would like to do a historic remark relating to esti- 
mation of the NPBF processes since in some works (e.g., see [36,35]) was 
expressed a surprise why these processes were not on a market during many 
years. First work [33], although found correct analytic expression for ISq 
pairing of neutrons, numerically underestimated the cmissivity by an order 
of magnitude. Also there was no statement on the dominance of the pro- 
cess in the cooling scenario (i.e. over the MU). Asymptotic behaviour of 
the emissivity e[nPBF] - 10 20 T 9 7 exp(-2zi/r) for T < A, T 9 = T/10 9 K, 
as follows from expression (lb) of [33] and from their rough asymptotic 
estimate of the integral (see below (13b)), shows up nor a large one nu- 
cleon phase space factor (<~ 10 28 ) nor appropriate temperature behaviour. 
Namely this underestimation of the rate in [33] (we again point out that 
analytic expression (la) is correct) and absence of mentioning of possible 
dominance of the process over MU, became the reason that this important 
result was overlooked during many years. Reference [24] overlooked sign of 
anomalous diagram for S pairing and included ex g 2 A term (g^ is axial-vector 
coupling constant) which should absent for S pairing giving, nevertheless, 
the main contribution for the 3P2 pairing case. A reasonable numerical es- 
timate of the emissivity was presented valid for both S and P pairings in- 
cluding NN correlation effects into consideration. Uncertainty of this esti- 
mate is given by a factor (0.5 -J- 2) which is allowed by variation of not too 
well known correlation factors. Namely this estimate was then incorporated 
within the cooling code in [26]. Correct asymptotic behaviour of the emis- 
sivity is e[nPFB,pPFB] ~ 10 28 (Z\/McV) 7 (T//i) 1 / 2 cxp(-2Z\/r) for T < A, 
that showed up a huge one-nucleon phase space factor and very moderate 
T dependence of the pre-factor. Thereby the value of the rate was related 
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in [24] to the value of the pairing gap. The possibility of the dominant role 
of this process (even compared to enhanced MMU and PU rather than only 
MU) was unambiguously stressed. Unfortunately [24] had a number of obvi- 
ous misprints which were partially corrected in subsequent papers. Although 
it does not excuse the authors these misprints can be easily treated by an 
attentive reader. Importance of NPFB processes was then once more stressed 
in review [16]. Reference [26] was the first that quoted the previous result [33]. 
It incorporated most important in-medium effects in the cooling code, among 
them nPBF and pPBF as equally important processes. However importance 
of pPBF process was then overlooked in the subsequent papers where emis- 
sivity of this process was several times incorrectly reproduced. Work [36] 
supported conclusion of [26] on importance of NPBF processes as governing 
the cooling scenario at T < T c . 

The medium modifications of all the above mentioned rates result in a 
pronounced density dependence (for NPBF processes mainly via dependence 
of the pairing gaps on the density and dependence on the NN correlation 
factors), which links the cooling behavior of a neutron star decisively to its 
mass [22,23,16,26]. As the result, the above mentioned medium modifications 
lead to a more rapid cooling than obtained in the "standard scenario ". Hence 
they provide a possible explanation for the observed deviations of some of the 
pulsar temperatures from the "standard" cooling. Particularly, they provide 
a smooth transition from "standard" to "non-standard" cooling for increasing 
central star densities, i.e., star masses. Thus by means of taking into account 
of most important in-medium effects in the reaction rates one is indeed able 
to achieve an appropriate agreement with both the high as well as the low 
observed pulsar temperatures that leads to the new "nuclear medium cooling 
scenario". Using a collection of modern EQS for nuclear matter, which cov- 
ered both relativistic as well as non-relativistic models, [26] also has demon- 
strated a relative robustness of these in-medium cooling mechanisms against 
variations in the EQS of dense NS matter. 

At initial stage (T > T opac ) the newly formed hot NS is opaque for neu- 
trinos/antineutrinos. Within FOPE model the value T opac was estimated in 
[6]. Elastic scatterings were included in [37,38] and pion condensation effect 
on the opacity was discussed in [39] . Medium effects dramatically affect the 
neutrino/ antineutrino mean free paths, since \ v (p) oc 1/ | M | 2 , where M 
is the reaction matrix element. Thereby, one-nucleon elastic scattering pro- 
cesses, like Nv — ► Nv 1 for energy and momentum transfer u> < qvFN are 
suppressed by NN correlations [24,16,40,41]. Neutrino/ antineutrino absorp- 
tion in two-nucleon MMU and MNB processes is substantially increased with 
the density (since | M | 2 for MMU and MNB processes increases with the 
density) [23,24,16]. Thus more massive NS are opaque for neutrinos up to 
lower temperatures that also results in a delay of neutrino pulse. Within the 
QPA for the nucleons the value T opac was estimated with taking into ac- 
count of medium effects in [23,16]. References [42,43,16] considered possible 
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consequences of such a delay for supernova explosions. On the other hand, 
at T > (1 -r- 2)MeV one should take care of the neutrino/antineutrino ra- 
diation in multiple NN scatterings (Landau Pomcranchuk Migdal (LPM) 
effect) when averaged neutrino-antineutrino energy, u> v a ~ several T, be- 
comes to be smaller than the nucleon width TV [28]. Numerical evaluations 
of -T/v in application to MNB processes were done in [44] and [45] within the 
Bruckner scheme and the Bethe-Salpeter equation, respectively. The LPM 
effect suppresses the rates of the neutrino elastic scattering processes on nu- 
cleons and also it suppresses MNB rates. For NS of rather low mass ( Ss Mq ) 
the suppression of the rates of neutral current processes due to the multiple 
collision coherence effect prevails over the enhancement due to the pion soft- 
ening, and for sufficiently massive NS (<; 1.4Mq) the enhancement prevails 
the suppression. MMU cmissivity remains to be almost unaffected by the 
LPM effect since averaged De energy ~ pp e is rather large (3> Jjv). 

The paper is organized as follows. Sect. 2 discusses basic ideas of the 
Fermi liquid approach to description of nuclear matter. The NN interaction 
amplitude is constructed with an explicit treatment of long-ranged soft pion 
mode and vertex renormalizations due to NN correlations. The meaning of 
the pion softening effect is clarified and a comparison of MOPE and FOPE 
models is given. Also renormalization of the weak interaction in NS matter is 
performed. Sect. 3 discusses the cooling of NS at T < T opac . Comparison of 
cmissivities of MMU and MU processes shows a significant enhancement of in- 
medium rates. Then we discuss DU-like processes and demonstrate medium 
effect due to vertex renormalizations. The role of in-medium mechanisms in 
the cooling evolution of NS is then demonstrated within a realistic cooling 
code. Next we consider influence of in-medium effects on the neutrino mean 
free path at initial stage of NS cooling. Essential role of multiple NN colli- 
sions is discussed. Sect. 4 presents OTF in non-equilibrium closed diagram 
technique in the framework of QPA for the nucleons and also beyond the 
QPA incorporating genuine particle width effects. 

2 Nuclear Fermi liquid description 

2.1 NN interaction. Separation of hard and soft modes 

At temperatures of our interest (T <C £fk) neutrons are only slightly excited 
above their Fermi sea and all the processes occur in a narrow vicinity of 
£Fn- In such a situation Fermi liquid approach seems to be the most efficient 
one. Within this approach the long-ranged diagrams are treated explicitly 
whereas short-scale diagrams are supposed to be the local quantities given 
by phenomenological so called Landau-Migdal (LM) parameters. Thus using 
argumentation of Fermi liquid theory [46-48,16] the retarded JVJV interaction 
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amplitude is presented as follows (see also [49]) 

X=X+XX+>CX> 

where 

(2) 



x = x + y^\ 



The solid line presents the nucleon, whereas double-line, the A isobar. The 
double- wavy line corresponds to the exchange of the free pion with inclusion 
of the contributions of the residual S wave ttNN interaction and 7T7t scat- 
tering, i.e. the residual irreducible interaction to the nucleon particle-holes 
and delta-nucleon holes. The full particle-hole, delta-nucleon hole and pion 
irreducible block (first block in (2)) is by its construction essentially more 
local then contributions given by explicitly presented graphs. Thereby, it is 
parameterized with the help of the LM parameters. In the standard Landau 
Fermi liquid theory fermions are supposed to be at their Fermi surface and 
the Landau parameters are further expanded in Lcgcndrc polynomials in the 
angle between fermionic momenta. Luckily, only zero and first harmonics en- 
ter physical quantities. The momentum dependence of the residual part of 
nuclear forces is expected to be not as pronounced and one can avoid this 
expansion. Then these parameters, i.e. /„„, f np and g nn , g n p in scalar and 
spin channels respectively, are considered as slightly momentum dependent 
quantities. In principle, they should be calculated as functions of the density, 
neutron and proton concentrations, energy and momentum but, simplifying, 
one can extract them from analysis of experimental data. 

The part of interaction involving A isobar is analogously constructed 



(3) 



The main part of the N A interaction is due to the pion exchange. Although 
information on local part of the N A interaction is rather scarce one can 
conclude [16,19] that the corresponding LM parameters are essentially smaller 
then those for NN interaction. Besides, at small transferred energies u) <C m^ 
the Z\-nucleon hole contribution is a smooth function of to and k in difference 
with the nucleon-nucleon hole (NN -1 ) contribution. Therefore and also for 
simplicity we neglect the first graph in r.h.s. of (3). 

Straightforward resummation of (1) in neutral channel yields [24,16] 


C - jf = Co n + Z%a x ■ <r 2 ) + fl N T^{a x ■ k)(a 2 ■ k), (4) 
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?% = faplifafl), Znn = 9nnl(.9nn), %np = 9npl(.9nn), a,0={n,p), 
%?n = J 2 (9nn)D%,, T^p = - 7 pp 7( £/„„)£>$, , T^ = ^ pp D^o, 

1- 1 (x) = l-2xC A* n , lpp = (1 - 4 5 Co4? n h(ff7m), (5) 

/nn = fpp = f + f, fnp = f - f, 9nn = #pp = .9 + #', and ,g np = g - g', di- 
mensional normalization factor is usually taken to be Co = tt 2 /[mpfPF(go)] — 
300 MeV- fm 3 ~ 0.77m- 2 , D* is the full retarded Green function of tt°, A af3 
is the corresponding NN^ 1 loop (without spin degeneracy factor 2) 



A ^ = O ' Ann{uj ~ q) ~ m « 2(47r2)_1 ( ln \ ± ^ ~ 2vFn ) ' (6) 



A nn ~ — m*pp' n (27r 2 )~ 1 , for w <C qvp n ,q <C 2pFn, and we for simplicity 
neglect proton hole contributions due to a small concentration of protons. 
Resummation of (1) in the charged channel yields 



X 



C P = ■ =Ca[Tg + zZ p( T 1 - ( T 2 )+f 2 N T n R p (<T 1 -k)(<T 2 -k), (7) 



V 



T* p = 2/'7(/'), -2& = 2.9 W), ^ = 7V)££- , (8) 

j- 1 (x) = l-4xC A^ p . 

The LM parameters are rather unknown for isospin asymmetric nuclear mat- 
ter and for g > qq. Although some evaluations of these quantities have been 
done, much work is still needed to get convincing results. Therefore for es- 
timates wc will use the values extracted from atomic nucleus experiments. 
Using argumentation of a relative locality of these quantities we will suppose 
the LM parameters to be independent on the density for g > Qq. One then 
can expect that the most uncertain will be the value of the scalar constant 
/ due to essential role of the medium-heavy a meson in this channel. But 
this parameter does not enter the tensor force channel being most important 
in our case. Unfortunately, there are also essential uncertainties in numerical 
values of some of the LM parameters even for atomic nuclei. These uncertain- 
ties are, mainly, due to attempts to get the best fit of experimental data in 
each concrete case slightly modifying parameterization used for the residual 
part of the NN interaction. E.g., calculations [47] gave / ~ 0.25, /' ~ 1, 
g ~ 0.5, g' ~ 1 whereas [50-52], including QP renormalization pre-factors, 
derived the values / ~ 0, /' ~ 0.5 + 0.6, g ~ 0.05 ± 0.1, g' ~ 1.1 ± 0.1. 

Typical energies and momenta entering NN interaction of our interest 
are to ~ and k ~ PFn- Then rough estimation yields ^(gnntW — 0, fc ~ 
PFn, Q = Qo) — 0.35 -j- 0.45. For u) — k ~ T typical for the weak processes 
with participation of vv one has r )~ l (g n n, ui — k ~ T, g — go) — 0.8 -j- 0.9. 
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2.2 Virtual Pion Mode 

Resummation of diagrams yields the following Dyson equation for pions 



■v^/ ^^<\» + -v^v^ J^>\, 



JJH 

res 



(9) 
The ttNA full-dot-vertex includes a background correction due to the pres- 
ence of the higher resonances, II^ es is the residual retarded pion self-energy 
that includes the contribution of all the diagrams which are not presented ex- 
plicitly in (9), as S wave ttNN and 7T7t scatterings (included by double- wavy 
line in (2)). The full vertex takes into account NN correlations 



^►~ 





(10) 



Due to that the nucleon particle- hole part of il^o is ex "f(g n n) and the nucleon 
particle-hole part of II 7r ± is ex j(g'). The value of the TV TV interaction in the 
pion channel is determined by the full pion propagator at small u> and k ~ 
PF n , i.e. by the quantity uj 2 (k) = —{D^)~ 1 {uj = 0,k,/j,„). Typical momenta 
of our interest are k ~ PFn- Indeed the momenta entering the NN interaction 
in MU and MMU processes are k — pp n i the momenta governing the MNB 
are k = kg [23] where the value k = fco ~ (0.9 H- l)pFn corresponds to the 
minimum of uj 2 (k). The quantity u) = w(fco) has the meaning of the effective 
pion gap. It is different for ir° and for ^ since neutral and charged channels 
are characterized by different diagrams permitted by charge conservation, 
thus also depending on the value of the pion chemical potential, fi^+ =/= /i v - ^ 
0, /i„-o = 0. For T <C e_Fn,e_F P , H„- — He = £_Fn — £Fp, as follows from 
equilibrium conditions for the reactions n — * pn~ and n — * pev. 

Change the sign of lo 2 symbolizes the pion condensate phase transition. 
Typical density behaviour of u> 2 is shown in Fig. 3. At g < (0.5 -j- 0.7)^oj °nc 
has u 2 = m 2 , — /j, 2 . For such densities the value uj 2 (pF n ) essentially deviates 
from m 2 , — /zj tending to m 2 + p 2 Fn — [i 2 in small density limit. 

At the critical point of the pion condensation (g = g C7r ) the value uj 2 with 
artificially switched off irir fluctuations (dashed line in Fig. 3) changes its sign. 
In reality -kit fluctuations are significant in the vicinity of the critical point 
and there occurs the first-order phase transition to the inhomogeneous pion- 
condensate state [53-55]. Thereby there are two branches (solid curves in Fig. 
3) with positive and respectively negative values for uj 2 . Calculations of [55] 
demonstrated that at g > g c1T the free energy of the state with ui 2 > 0, where 
the pion mean field is zero, becomes larger than that of the corresponding 
state with lo 2 < and a finite mean field. Therefore at g > g C7r the state 
with uj 2 > is metastable and the state with lo 2 < and the pion mean field 
(fin =/= becomes the ground state. 

The quantity u> 2 demonstrates how much the virtual (particle-hole) mode 
with pion quantum numbers is softened at given density. The ratio a = 
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Fig. 3. Effective pion gap (for (i^a = 0) versus baryon density, see [16]. 

D 7r [med.]/-D 7 r[vac.] ~ 6 for g = go, lo = 0, k = pfn and for isospin symmetric 
nuclear matter. However this essential so called "pion softening" [48] does not 
significantly enhance the NN scattering cross section due to a simultaneous 
essential suppression of the irNN vertex by NN correlations. Indeed, the 
ratio of the NN cross sections calculated with FOPE and MOPE is 



R 



a[FOPE] 7 V, to ~ 0, k ~ PFN )(ml + p 2 FN f 



(11) 



a [MOPE] " to 4 <J>fn) 

and for g = go we have R ^ 1, whereas for g — 2 go we already get R ~ 10. 

As follows from numerical estimates of different 7 factors entering (4) and 
(7), the main contribution to NN interaction for g > go is given by MOPE 



X 



(12) 



whether this channel (T ex {<j\ ■ fc)(er 2 • fe)) of the reaction is not forbidden or 
suppressed by some specific reasons like symmetry, small momentum transfer, 
etc. The g meson contribution to NN interaction is partially included in g a p, 
other part contributing to T and T is minor (oc uj 2 /m e 2 ). Thus instead of 
FOPE+0 exchange, as the model of NN interaction which has been used in 
[6] in calculation of the cmissivities of the two-nucleon reactions within the 
"standard scenario" of NS cooling one should use the full NN interaction 
given by (4) and (7) or, simplifying, approximated by its MOPE part. 



2.3 Renormalization of the weak interaction. 

The full weak coupling vertex that takes into account NN correlations is 
determined by (10) where now the wavy line should be replaced by the lepton 
pair. Thus for the vertex of our interest, Ni —* N2W, we obtain [24,16] 

G 



V, 



P 



V2 



[7(/')*o - 9AW)b 



(13) 



12 D.N. Voskresensky 

for the (3 decay and 
G 



V 71.71, 



2V2 



[l{fnn)k ~ 9Al{9nn)l<- 



V N = — 

PP 2^2 



[k pp 1 - gAjpplcr] j(14) 



Kpp = Cv — 1fnpl{fnn)CoA nn , j pp = (1 - 4gC A nn ) j(g nn ), (15) 

for processes on the neutral currents NiN 2 — ► N1N2VV, V pp = VJl + VS 1 , G ~ 



1.17- 10" 5 GeV 

,2 



2 is the Fermi weak coupling constant, cy = 1 — 4 sin #w> 
iy ~ 0.23, qa — 1-26 is the axial-vector coupling constant, and l^ = 



sm 

u(gi)7 (U (l— 75)u((Z2) is the lepton current. The tensor force contribution ~ q 2 
is small for typical | q |~ T or pFe, and for simplicity is omitted. Possibility 
of Brown-Rho scaling can be easily incorporated by scaling of qa [56] . With a 
decrease of g A the value g C7T increases. But it remains finite (~ 2go according 
to [57]) due to attractive contribution of the A isobar in (9), whereas one 
would expect g C7r — ► oo for g A — > 1 ignoring A contribution. 

The 7 factors renormalize the corresponding vacuum vertices. These fac- 
tors are essentially different for different processes involved. The matrix el- 
ements of the neutrino/ antineutrino scattering processes Nv — > Nv and of 
MNB behave differently in dependence on the energy-momentum transfer 
and whether N = n or N = p in the weak coupling vertex. Vertices 




-< 



(16) 



N N N 

are modified by the correlation factors (5) and (8). For N — n these are 
7(grm, w, q) and "/(f nn ,u;,q) leading to an enhancement of the cross sections 
for lu > qvFn and to a suppression for lu < qvp n . Rcnormalization of the 
proton vertex (vector part of V p ^ + V^ p ) is governed by the processes [24,58] 



v + ...(17) 



being forbidden in vacuum. For the systems with lSo proton-proton pairing, 
oc g A contribution to the squared matrix element (see (14)) is compensated by 
the corresponding contribution of the diagram with anomalous Green func- 
tions of protons. The vector current term is oc Cy in vacuum whereas it is 
ex K pp in medium (by first graph (17)). Thereby the corresponding vertices 
with participation of proton are enhanced in medium compared to the small 





vacuum value (ex 



0.006) leading to enhancement of the cross sections 
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(up to ~ 10 4- 10 2 times for 1.5 ~ 3po depending on parameter choice). Also 
enhancement factor (up to ~ 10 2 ) comes from the intermediate in-medium 
photon (j m ) whose propagator contains 1/to 2 ~ 1/e 2 , where m 7 is the ef- 
fective spectrum gap, that compensates small e 2 factor from electromagnetic 
vertices [58]. It is included by replacement V^ — > V pp . Other processes per- 
mitted in intermediate states like processes with pp~ x and with the pion are 
suppressed by a small proton density and by q 2 ~ T 2 pre-factors, respectively. 
First diagram (17) was considered in [24] where the pPBF process was sug- 
gested, and then in [16,26], and it was shown that nPBF and pPBF processes 
may give contributions of the same order of magnitude. Several subsequent 
papers overlooked these results and rediscovered the pPBF process ignor- 
ing the nucleon and the electron correlation effects. Specific contribution of 
second diagram for pPFB process was recently discussed in [59] . 

Paper [27] gives an other example demonstrating that although the vac- 
uum branching ratio of a kaon decays is r{K~ — > e~ + i> e )/ T{K~ —* 
\i~ + Vfj) s=a 2.5 x 10~ 5 in medium (due to Ap^ 1 decays of virtual K~) it 
becomes to be of order of unit. Thus we see that in dependence of what re- 
action channel is considered in-medium effects may or strongly enhance the 
reaction rate under consideration or substantially suppress it. 



2.4 Inconsistencies of FOPE model 

Since FOPE model became the base of the "standard scenario" for cooling 
simulations we would like first to demonstrate principal inconsistencies of the 
model for the description of interactions in dense (g > go) baryon medium. 
The only diagram in FOPE model which contributes to the MU and NB is 




(18) 



firN 

Dots symbolize FOPE. This is first available Born approximation diagram, 
i.e. second order perturbative contribution in f^N- In order to be theoretically 
consistent one should use perturbation theory up to the very same second 
order in f n ^ f° r a U the quantities. E.g., pion spectrum is then determined 
by pion polarization operator expanded up to the very same order in f„jy 

fnN, 

u? ~ m 2 +k 2 +n* V, k, g), iT fl (w, k,g)= ■ ■ 

%N ( 19 ) 

The value 12 (u>, k, g) is easily calculated containing no any uncertain param- 
eters. For u — > 0, k ~ pf of our interest and for isospin symmetric matter 



o 



II« a - ao - i/Joc, Qo c 2m >™?Z» > 0, ft a UMIIn > o. (20) 
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Replacing this value to (19) we obtain a solution with \u> < already for 
g > 0.3go that would mean appearance of the pion condensation. Indeed, 
the mean field begins to increase with the time passage <p ~ exp(— \wt) ~ 
exp(crf//3) until repulsive irir interaction will not stop its growth. But it is 
experimentally proven that there is no pion condensation in atomic nuclei, i.e. 
even at g = go- The puzzle is solved as follows. FOPE model does not work 
for such densities. One should replace FOPE by the full NN interaction given 
by (4) , (7) . Essential part of this interaction is due to MOPE with vertices 
corrected by NN correlations. Also the NN -1 part of the pion polarization 
operator is corrected by NN correlations. Thus 



^-f~^^- - n o{u,k,Q)7{g',ui,k,Q) 



(21) 



being suppressed by the factor j(g',u — 0, k ~ pf, Q — go) — 0.35 -j- 0.45. 
Final solution of the dispersion relation (19), now with full II instead of 77°, 
yields uo > for g = go whereas the solution with iu> < 0, which shows the 
begining of pion condensation, appears only for g > g C7r > go- 

3 Neutrino cooling of neutron stars 

3.1 Emissivity of MMU process 

Since DU process is forbidden up to sufficiently high density Qaj, the main 
contribution for g < g c \j and T opac > T > T c comes from MMU pro- 
cesses schematically presented by two diagrams of Fig. 1. MNB reactions 
give smaller contribution [23] . For densities g <C go the main part of the NN 
interaction amplitude is given by the residual NN interaction. In this case 
the NN interaction amplitude can be better treated within the T matrix 
approach which sums up the ladder diagrams in the particle-particle chan- 
nel rather then by LM parameters. Calculations of MNB processes with the 
vacuum T matrix [60] found essentially smaller emissivity then that given 
by FOPE. Also the in-medium modifications of the T matrix additionally 
suppress the rates of both MMU and MNB processes, see [61]. Thus even at 
small densities the FOPE model may give only a rough estimate of the emis- 
sivity of two nucleon processes. At g £ (0.5 -j- 0.7) go reactions in particle-hole 
channel and more specifically with participation of the soft pion mode begin 
to dominate. 

Evaluations [23,42,43,16] showed that the dominating contribution to 
MMU rate comes from the second diagram of Fig. 1, namely from contri- 
butions to it given by the first two diagrams of the series 
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e 




n n n n (22) 

whereas the third diagram, which naturally generalizes the corresponding 
MU(FOPE) contribution, gives only a small correction for g > go. The emis- 
sivity from the two first diagrams in a simplified notation [16,26] reads 



MMU 



MOPE] 



2.4 • 10 24 T$ ( ± 



1 4 



aCJ w ±(pFn) 



r s F 1 Q{A r , 



.Qo 

C(A> 



10/3 



(m* 



''N 



1 4 



aQ v o (pFn) 



erg 



(23) 



where Tg = T/10 9 K is the temperature, m* and m* are the nonrelativistic 
effective neutron and proton masses, and the correlation factor r s is roughly 

* ~' 7^(^ - PFe,q - PFeh 2 (gnn,U ~ 0, fc = p Fn )j 4 (g' ,0, p Fn ) , 

7 2 (/>,9) + 3.9i7 2 (5>,9) 



r 



7 /3(w,g) 



l + 3 5 i 
and the second term in the factor 



F, 



1 



PFe,q — PFe) 



2/3 



(24) 



(25) 



4 7 2 ( 5 ',0,pFn)7|(^ 

is the contribution of the pion decay from intermediate states (first diagram 
(22)). The quantity r effectively accounts for an appropriate product of the 
NN correlation factors in different ttNiN 2 vertices. For charged pions the 
value [in ^ is incorporated in the expression for the effective pion gap, for 
neutral pions fi n — 0. The value a ~ 1 depends on condensate structure, 
a = 1 for q < g C7r , and a = v2 taking account of the new excitations on the 
ground of the charged 7r condensate vacuum for g > g C7! . The factor 

t(A \ ~ / exp (~ A n/T) T<T cN , . . 

« An ^\i T>T cN , N={n,p) ^ 

estimates the suppression caused by the nn and pp pairings. Deviation of 
these factors from simple exponents can be incorporated as in [35]. 

The ratio of the cmissivities of MMU(MOPE) and MU(FOPE) is roughly 

£ MMU [M0PE] ^ 3 j2(g nn ,0,p Fn )^(gi,0,p Fn ) 



~MU 



[FOPE] 



10 J 



(o/eo) 



10/3 



(27) 



U!^o(p Fn )ujl ± (p Fn ) 

For g ~ go this ratio is ~ 10 whereas being estimated with the only third 
diagram (22) it would be less then unit. 
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3.2 Emissivity of DU-like processes 

NPBF processes. The one-nucleon processes with neutral currents given by 
the second diagram (16) for N — (n,p) are forbidden at T > T c by energy- 
momentum conservations but they can occur at T < T c . Then physically the 
processes relate to NPBF, see Fig. 2. However they need special techniques 
to be calculated [33,24]. These processes n —* nvv and p — > pvv play very 
important role in the cooling of superfluid NS, see [24,25,16,26,36,35]. The 
emissivity for 3 types of neutrinos is given by [24,25] 

f pBFl = 3-4G 2 (g l7 2 (/n„)+6ffA7 2 (ffn«))PFn<^ / A,\ 

e[n J 15tt 5 V T ) 

-C-10 28 ^) 1/3 ^(^) 7 /(^)^^, T<T cni (28) 
\QoJ m N \ MeV / \ T J cm d • sec 

where £i = 1, £2 = for S'-pairing and £1 = 2/3, £2 = 4/3 for P-pairing, 
compare [35]. I removed some misprints existed in [33,24,25,16]. For neu- 
trinos Lo(q) = q and correlations are not so essential as it would be for 
u> <C q. Taking j 2 ~ 1.3 in the range of 5*-pairing we get £ ~ 5 whereas 
for the P-pairing with g* A ~ 1.1 we obtain Q ~ 9, in agreement with nu- 
merical evaluations [24] (£ ~ 7) used then within the cooling code in [26], 
I(x) = r o °°ch 5 2/dy/(exp(xch2/) + l) 2 , I(x > 1) ~ cyLv{-2A/T)^TtT/AA, 
that serves an appropriate asymptotic (28). Emissivity of the process p — > pvv 
is given by [24] 

£(P PBF) = 12G2 «i«gp + ^j, + 6)p f pr n ;^ j (±) fT< Tcp[29) 



15tt 5 V T , 

and $2=0 for protons paired in S-state in NS matter, $3 ^ 1 is due to 
the second diagram (17) and has a complicated structure [58,59]. For the 
process (29) the part of TV TV and ee correlations is especially important. One 
has Kp p ~ 0.05 -7- 1 for 1 — 3qo and $3 <~ 1, and Kp p + $3^1 instead of a 
small Cy ~ 0.006 value in absence of correlations. Thereby, in agreement with 
[24,26,58,59], the emissivity of the process p — > pvv can be compatible with 
that for n — ► nvv in dependence on the relation between A p and A n . 

NPBF processes are very efficient for T <T C competing with MMU pro- 
cesses. The former win the content for not too massive stars. Analysis of 
above processes supports also our general conclusion on the crucial role of 
in-medium effects in the cooling scenario. 

Pion (kaon) condensate processes. The P wave pion condensate can 
be of three types: 7r+, 7r ± , and tt° with different values of the critical densities 
Qcir = (ft;7r± j C7r +i Qctt°)i see [48]. Thus above the threshold density for the 
pion condensation of the given type, the neutrino emissivity of the MMU 
process (23) is to be supplemented by the corresponding PU processes 
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n -*. 



z z z 

x e x e x iv 

x xx (30) 

The emissivity of the charged pion condensate processes with inclusion of the 
NN correlation effect (in a simplified treatment) renders, see [22,16], 

e[PU] ~ 7- 10 26 P£^!^<jl(p FeiPFe ) ^O,p Fn )lfsin 2 0-^I-(31) 
m^ m N ' cm d sec 



Here g > g C7r and sin6> ~ ^/2\ui 2 \/m 2 for 6 <C 1. Of the same order of 
magnitude are emissivities of other possible 7r condensate reactions, e.g. for 
7i7r" — > npev process at 9 <C 1 the numerical factor is about two times larger. 
Since 7r condensation probably reduces the energy gaps of the supcrfluid 
states by an order of magnitude, see [62], we may assume that superfluidity 
vanishes above g C7r . Finally we note that though the PU processes have gen- 
uinely one-nucleon phase-space volumes, their contribution to the emissivity 
is suppressed relative to the DU by an additional 7 2 (g' ,0,pFn) suppression 
factor due to existence of the extra (ttNN) vertex in the former case. 

Fig. 4 compares the mass dependence of the neutrino cooling rates L v /Cy, 
L v is the neutrino luminosity, CV is the heat capacity, for MMU(MOPE) 
and MU(FOPE) for non-superfluid matter. For the solid curves, the neutrino 
emissivity in pion-condensed matter is taken into account according to (31) 
and (23) with the parameter a = \/2. To be conservative we took g C7T ~ 
3go- Dashed curves correspond to the model where no pion condensation is 
allowed. As one sees, the medium polarization effects included in MMU may 
result in three order of magnitude increase of the cooling rate for the most 
massive stars. Even for stars of a rather low mass the cooling rate of MMU 
is still several times larger than for MU because even in this case the more 
efficient rate is given by the reactions shown by the right diagram in Fig. 1 
(first two diagrams of (22)). The cooling rates for the NS of M = 1.8M with 
and without pion condensate differ only moderately (by factor of 5 in this 
model). If we would take g cv to be smaller the ratio of emissivity PU to MMU 
would decrease and could even become 5s 1 in the vicinity of g CK . Opposite, 
the reaction rates for the FOPE model are rather independent of the star's 
mass for the stars with masses below the critical value 1.63Mq, at which 
transition into the pion condensed phase occurs, and then jump to the typical 
PU value. It is to be stressed that contrary to FOPE model, the MOPE model 
[23] consistently takes into account the pion softening effects for g < g CT: and 
both the pion condensation and pion softening effects on the ground of the 
condensate for g > g C7T . For g > g C K the kaon condensate processes come into 
play. Most popular is the idea of the S wave K~ condensation (e.g. see [13]) 
which is allowed at /i e > rn* K _ due to possibility of the reaction e — ► K~v. 
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• FOPE with Tt-cond. 

■ MOPE with ji-cond. 

■ MOPE without 71-cond 

• FOPE without 71-cond. 




Fig. 4. Cooling rate due to neutrino emission as a function of star mass for a 
representative temperature of T = 3 x 10 s K. Superfluidity is neglected. 



Analogous condition for the pions, fi e > m* _ , is not fulfilled owing to a strong 
S wave ttNN repulsion [48,16] (again in-medium effect!) otherwise S wave 
7r~ condensation would occur at smaller densities then K~ condensation. 
The neutrino emissivity of the K~ condensate processes is given by equation 
analogous to (31) with a different NN correlation factor and an additional 
suppression factor due to a small contribution of the Cabibbo angle. However 
qualitatively scenario that permits kaon condensate processes is analogous to 
that with the pion condensate processes. 

Other resonance processes. There are many other reaction channels 
allowed in the medium. E.g, any Fermi liquid permits propagation of zero 
sound excitations of different symmetry related to the pion and the quanta 
of a more local interaction determined here via / a>/ g and g a ,p- These excita- 
tions being present at T ^ may also participate in the neutrino reactions. 
The most essential contribution comes from the neutral current processes [23] 
given by first two diagrams of the series - 









(32) 



Here the dotted line is zero sound quantum of appropriate symmetry. These 
are the resonance processes (second, of DU-type) analogous to those processes 
going on the condensates with the only difference that the rates of reactions 
with zero sounds are proportional to the thermal occupations of the corre- 
sponding spectrum branches whereas the rates of the condensate processes 
are proportional to the modulus squared of condensate mean field. Contribu- 
tion of the resonance reactions is as a role rather small due to a small phase 
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space volume (q ~ T) associated with zero sounds. Please also bear in mind 
an analogy of the processes (32) with the corresponding phonon processes in 
the crust. 

DU processes. The proper DU processes in matter, as n — > pe~v e and 
per — ► nv e , 



'> d + ; (33) 



p p n 

should also be treated with the full vertices. They are forbidden up to the 
density g c u when triangle inequality pp n < pp p + pp e begins to fulfill. For 
traditional EQS like that given by the variational theory [18] DU processes 
are permitted for g > 5g»o- The emissivity of the DU processes renders 



z DU ~ 1.2 • 10 27 ^^ 



( loc^ev) ^ min K(AO,C(A0] n ^< (34) 



where /U; = /i e = /i M is the chemical potential of the leptons in MeV. In 
addition to the usually exploited result [15], (34) contains 7? pre-factor (13) 
due to NN correlations in the (3 decay vertices, see [24,26]. 

It was realized in [63] that the softening of the pion mode in dense neutron 
matter could also give rise to a rearrangement of singlc-fcrmion degrees of 
freedom due to violation of Pomcranchuk stability condition for g = g c p < 
g> C7r . It may result in organization of an extra Fermi sea for g c p < g < g C7T 
and at small momenta p < Q.2pp n , that in its turn opens a DU channel 
of neutrino cooling of NS from the corresponding layer. Due to the new 
feature of a temperature-dependent neutron effective mass, to* oc 1/T, we 
may anticipate extra essential enhancement of the DU process, corresponding 
to a reduction in the power of the temperature dependence from T 6 to T 5 . 
At early hot stage this layer becomes to be opaque for the neutrinos slowing 
the transport from the massive NS core to the exterior. 

Basing on the Brown-Rho scaling idea [56], we argued in [64] for the 
charged g meson condensation at a relevant density (g cg ~ 3go if to* would 
decrease to that density up to m g /2). If happend, it would open DU reaction 
already for g < g cg and close it for g > g cg due to an essential modification 
of the nuclear asymmetry energy. 

3.3 Comparison with soft X ray data 

The heat transport within the crust establishes homogeneous density profile 
at times ^ (1 -j- 10)yr. After that time the subsequent cooling is determined 
by simple relation CyT = —L, where Cy = J^ i Cy.i and L = J^ i Li are sums 
of the partial contributions to the heat capacity (specific heat integrated over 
the volume) and the luminosity (emissivity integrated over the volume). 
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3.0 
log(T/yr) 



Fig. 5. Cooling of non-suporfluid NS models of different masses constructed for 
the HV EOS [26]. Two graphs refer to cooling via MU(FOPE) +PU (left) and 
MMU(MOPE)+PU (right). In both cases, pion condensation is taken into account 
for the solid curves at g > g c -w, i.e. for M > 1.6Mq (g C7T is chosen to be 3go). The 
dashed curves in the right graph refer to a somewhat larger value of uj 2 (pFn), with- 
out pion condensation. The observed luminosities are labeled by dots. Possibilities 
of Fermi sea rearrangement and g condensation are ignored. 



The nucleon pairing gaps arc rather purely known. Therefore one may 
vary them. The "standard" and "nonstandard" scenarios of the cooling of NS 
of several selected masses for suppressed gaps are demonstrated in the left 
panel of Fig. 5, [26]. Depending on the star mass, the resulting photon lumi- 
nosities are basically either too high or too low compared to those given by 
observations. Situation changes, if the MMU process (23) is included. Now, 
the cooling rates vary smoothly with the star mass (see right panel of Fig. 
5) such that the gap between standard and non-standard cooling scenarios is 
washed out. More quantitatively, by means of varying the NS mass between 
(1 -j- 1.6) M, q , one achieves an agreement with a large number of observed 
data points. This is true for a wide range of choices of the uj 2 parameteriza- 
tion, independently whether pion (kaon) condensation can occur or not. Two 
paramcterizations presented in Fig. 5 with pion condensation for g > 3go 
and without it differ only in the range which is covered by the cooling curves. 
The only point which does not agree with the cooling curves belongs to the 
hottest pulsar PSR 1951+32. Other three points which according to Fig. 5, 
right panel, are also not fitted by the curves can be easily fitted by slight 
changes of the model parameters. The high luminosity of PSR 1951+32 may 
be due to internal heating processes, cf. [10]. 
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Fig. 6. Cooling of NS with different masses constructed for the HV EOS [26]. The 
cooling processes are MU-VS86, PU (only solid curves), NPBF, PPBF, and DU 
(only in the right graph). The dashed curves refer to the M — 1.7 and 1.9Mq 
models without pion condensate. Possibilities of Fermi sea rearrangement and g 
condensation are ignored. 



We turn now to cooling simulations where the MU, NPBF, DU and PU 
take place simultaneously. Parameters of the pairing gaps are from Fig. 6 of 
[26]. Fig. 6 shows the cooling tracks of stars of different masses, computed for 
the HV EQS. Very efficient at T < T c become to be NPBF processes which 
compete with MMU processes. The former prevail for not too massive stars in 
agreement with estimation [24] . The DU process is taken into account in the 
right graph, whereas it is neglected in the left graph. The solid curves refer 
again to the Cj 2 parameterization with phase transition to pion condensate, 
the dashed curves to the one without phase transition (see Fig. 4). For masses 
in the range between 1.0 and 1.6 Mq, the cooling curves pass through most of 
the data points. We again recognize a photon luminosity drop by more than 
two orders of magnitude for the 1.7 M Q mass star with pion condensate, due 
to suppression of the pairing gaps in this case. This drop is even larger if the 
DU is taken into account (right graph) . This allows to account for the photon 
luminosity of PSR 1929+10. 

Thus, comparison with the observed luminosities shows that one gets 
quite good agreement between theory and observations if one includes into 
consideration all available in-medium effects assuming that the masses of the 
pulsars are different. We point out that the description of these effects is 
constructed in essentially the same manner for all the hadron systems as NS, 
atomic nuclei and heavy ion collisions, cf. [16,49]. 
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3.4 Neutrino opacity 

Important part of in-medium effects for description of neutrino transport at 
initial stage of NS cooling was discussed in [23,42,43,16], where correlation ef- 
fects, pion softening and pion condensation, the latter for g > g C7T , were taken 
into account. The neutrino/antineutrino mean free paths can be evaluated 
from the corresponding kinetic equations via their widths -Op) = — 2Im77^_^ , 

where II R is the retarded self-energy, or within the QPA for the nucleons they 
can be also estimated via the squared matrix elements of the corresponding 
reactions. Thereby the processes which the most efficiently contribute to the 
emissivity are at early times (for T > 1 MeV) responsible for the opacity. 

In the above "nuclear medium cooling scenario" at T > T c the most es- 
sential contribution was from MMU. Taking into account of NN correlations 
in the strong coupling vertices of two-nucleon processes like MMU and MNB 
suppresses the rates, whereas the softening of the pion propagator essentially 
enhances them. For rather massive NS MOPE wins the competition. The 
mean free path of neutrino/antineurtino in MMU processes is determined 
from the same diagrams (22) as the emissivity. Its calculation (see (23)) with 
the two first diagrams yields 
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(35) 

From the relation A„ ~ R follows evaluation of T opac . With the only first 
diagram we get a simple estimate 
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For averaged value of the density g ~ go corresponding to a medium-heavy 
NS (< 1.4M ) with Cj 2 {p Fn ) ~ O.Sml , 7 ~ 7 ~ (0.3 H- 0.4) we get T opac ~ 
(1 -j- 1.5) MeV that is smaller then the value T opac ~ 2 MeV estimated with 
FOPE [6]. For g ~ 2 go that corresponds to a more massive NS we evaluate 
T opac ~ (0.3 -j- 0.5) MeV. Thus pion softening results in a substantial decrease 
of neutrino/antineutrino mean free paths and the value of T opac . 

The diffusion equation determines the characteristic time scale for the 
heat transport of neutrinos from the hot zone to the star surface to ~ 
R 2 C\/o-~ 1 T~ 3 j \ v (a is the Stefan-Boltzmann constant), whereas it follows 
that t ~10 min. for T ~ 10 MeV and g ~ g , lo 2 ~ 0.8m,,, T ~ 0.4, 
m* N jm^s ~ 0.9, and to becomes as large as several hours for g ~ (2 -f- 3)go- 
These estimates demonstrate that more massive NS cool down more slowly 
at T > T opac and faster at subsequent times then the less massive stars. 

Due to in-medium effects neutrino scattering cross section on the neutrons 
shown by the first diagram (16) requires the NN correlation factor 

2/ x J 2 (fnn,u,q) + 3g 2 A -y 2 (g nn ,UJ,q) 
7„>, 9) YT ^ , (37) 
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as follows from (14). This results in a suppression of of the cross sections for 
lo < qvp n and in an enhancement for u > qvp n . Neutrino scattering cross 
sections on the protons are modified by 

2, x K 2 (f np J nn ,u;,q) + 3g 2 A rf(g nn ,uj,q) 

V(«.«) = YTW A ' (38) 

that results in the same order of magnitude correction as given by (37). 

Also there is a suppression of the vN scattering and MNB reaction rates 
for soft neutrinos (u> Ss (3 4- 6)T) due to multiple NN collisions 



I 



^z &» 



(LPM effect). One may estimate these effects simply multiplying squared 
matrix elements of the vN scattering and the MNB processes by the cor- 
responding suppression pre-factors [28]. Qualitatively one may use a general 
pre-factor Cq{lo) = uj 2 /[lj 2 +T^], where Jjv is the nucleon width and the lo is 
the energy of v or vv pair. In some works, see [65,66], correction factor, like 
Co, was suggested at an ansatz level. Actually one does not need any ansatze 
reductions. OTF, see [28], allows to calculate the rates using the exact sum 
rule. The modification of the charged current processes due to LPM effect is 
unimportant since the corresponding value of lo is ~ pf & S> -Tat. 

The main physical result we discussed is that in medium the reaction rates 
are essentially modified. A suppression arises due to NN correlations (for 
lo <C kvFN) and infra-red pre-factors (coherence effects), and an enhancement 
due to the pion softening (and NN correlations for w ~ q) and due to opening 
up of new efficient reaction channels. The pion softening demonstrates that 
already for densities g < go the nucleon system begins to feel that it may 
have 7r condensate phase transition for g > g c1T , although this g C7Z value might 
be essentially larger than go or even not achieved. 

4 The rate of radiation from dense medium. OTF 

Perturbative diagramss are obviously irrelevant for calculation of in-medium 
processes and one should deal with dressed Green functions. The QPA for 
fcrmions is applicable if the fermion width is much less than all the typical 
energy scales essential in the problem (/> <C u) c h)- In calculation of the emis- 
sivities of vv reactions the minimal scale is u c h — 6T, averaged vv energy for 
MNB reactions. For MMU u> c h — PFe- For radiation of soft quanta of fixed 
energy ui < T, lo c ] x ~ lo. Within the QPA for fcrmions, the reaction rate with 
participation of the fermion and the boson is given by [22,23] 
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For equilibrium (T^O) system there is the exact relation 



(< (plfo >)(p) = iD-++ | 



A, 



¥>c 



cxp(^) - 1 



'■Pc 



A B = -2ImD B \ 




(41) 

where (< 02<Pi >)(p) means the Fourier component of the corresponding 
non-equilibrium Green function and <p c is the mean field. Thus the rate of 
the reaction is related to the boson spectral function Ab and the width {Tb) 
being determined by the corresponding Dyson equation, see (9). Ab is the 
delta-function at the spectrum branches related to resonance processes, like 
zero sound. The poles associated with the upper branches do not contribute 
at small temperatures due to a tiny thermal population of those branches. 
There is also a contribution to ImD R proportional to ImII R given by the 
particle-hole diagram. Within the QPA taking Im part means the cut of the 
diagram. Thus we show [23] that this contribution is the same as that could 
be calculated with the help of the squared matrix element of the two-nucleon 
process ^^ 2 



(42) 



This is precisely what one could expect using optical theorem. Thus unlim- 
ited series of all possible diagrams with in-medium Green functions (see (16), 
(30), (32), (33)) together with two-fermion diagrams (as given by (22)) and 
multiplc-fcrmion diagrams (like (39)) would lead us to a double counting. The 
reason is that permitting the boson width effects (and beyond the QPA for 
fcrmions also permitting finite fcrmion widths) the difference between one- 
fermion, two-fermion and multiplc-fcrmion processes in medium is smeared 
out. All the states are allowed and participate in production and absorption 
processes. Staying with the QP picture for fermions, the easiest way to avoid 
mentioned double counting is to calculate the reaction rates according to 
(40), i.e. with the help of the diagrams of the DU-like type, which already 
include all the contributions of the two-nucleon origin. Multiple NN collision 
processes should be added separately. On the other hand, it is rather inconve- 
nient to explicitly treat all one-nucleon processes dealing with different spe- 
cific quanta instead of using of the full NN interaction amplitude whenever 
it is possible. Besides, as we have mentioned, consideration of open fermion 
legs is only possible within the QPA for fcrmions since Feynman technique 
is not applicable if Green functions of ingoing and outgoing fcrmions have 
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widths (that in another language means possibility of additional processes) . 
Thus the idea came [24,28] to integrate over all in-medium states allowing all 
possible processes instead of specifying different special reaction channels. 

In [24,28,27] it was shown that OTF in terms of full non-equilibrium 
Green functions is an efficient tool to calculate the reaction rates including 
finite particle widths and other in-medium effects. Applying this approach, 
e.g., to the antineutrino-lepton (electron, pT meson, or neutrino) production 
[24] we can express the transition probability in a direct reaction in terms of 
the evolution operator S, 

;tot 



JWxlvi _ (l - m)dgfd4 

dt (27r) 6 4w;wp 



^<0\S^\vl + X><Dl + X\S\0>, (43) 

{*} 



where we presented explicitly the phase-space volume of Pi states; lepton oc- 
cupations of given spin, ni , are put zero for v and v which are supposed to be 
radiated directly from the system (for T < T opac ). The bar denotes statistical 
averaging. The summation goes over complete set of all possible intermediate 
states {X} constrained by the energy-momentum conservation. It was also 
supposed that electrons/muons can be treated in the QPA, i.e. with zero 
widths. Then there is no need (although possible) to consider them in inter- 
mediate reaction states. Making use of the smallness of the weak coupling, 
we expand the evolution operator as S w 1 — i J_ T {Vw(x) S nuc \(x)}dxo , 
where Vw is the Hamiltonian of the weak interaction in the interaction repre- 
sentation, S nuc i is the part of the S matrix corresponding to the nuclear inter- 
action, and T{...} is the chronological ordering operator. After substitution 
into (43) and averaging over the arbitrary non-equilibrium state of a nuclear 
system, there appear chronologically ordered (G ), anti-chronologically or- 
dered (G ++ ) and disordered (G H and G~ + ) exact Green functions. 

In graphical form the general expression for the probability of the lepton 
(electron, muon, neutrino) and anti-neutrino production is as follows 

I v 




I 

representing the sum of all closed diagrams (—ill h ) containing at least one 
( — h) exact Green function. The latter quantity is especially important. Var- 
ious contributions from {X} can be classified according to the number N of 
G h lines in the diagram 

v I 
dWlf = (1 - ntidPqpdqf 
dt (27r) 6 4uj i? 0Ji 

This procedure suggested in [24] is actually very helpful especially if the 
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QPA holds for the fermions. Then contributions of specific processes con- 
tained in a closed diagram can be made visible by cutting the diagrams 
over the (H — ), ( — h) lines. In the framework of the QPA for the fermions 
G h = 2i:\np8{e + u — e° — HjeS R (e + n, p)) (np are fermionic occupations, 
for equilibrium np — l/[exp((e — /j,p)/T) + 1]), and the cut eliminating the 
energy integral thus requires clear physical meaning. This way one establishes 
the correspondence between closed diagrams and usual Feynman amplitudes 
although in general case of finite fcrmion width the cut has only a sym- 
bolic meaning. Next advantage is that in the QPA any extra G '", since it 
is proportional to np, brings a small (T/ep) 2 factor to the emissivity of the 
process. Dealing with small temperatures one can restrict by the diagrams of 
the lowest order in (G h G^ ), not forbidden by energy- momentum conser- 
vations, putting T = in all G ++ and G~~ Green functions. Each diagram 
in (44) represents a whole class of perturbative diagrams of any order in the 
interaction strength and in the number of loops. 

Proceeding further we may explicitly decompose the first term in (44) as 
I v I v I v 




XX XX 



(45) 



I v I 

The full vertex in the diagram (45) of given sign is irreducible with respect to 

the (-1 — ) and ( — h) nucleon-nucleon hole lines. This means it contains only 

the lines of given sign, all ( ) or (H — \-). Second diagram with anomalous 

Green functions exists only for systems with pairing. In the framework of 
the QPA namely these diagrams determine the proper DU and also NPBF 
processes calculated in [24,25] within OTF. In the QP picture the contribution 
to DU process vanishes for g < g c jj. Then the second term (44) comes into 
play which within the same QP picture contains two-nuclcon processes with 
one (G h G H ) loop in intermediate states, etc. 

The full set of diagrams for 77 *" can be further explicitly decomposed as 
series [28] (from now in brief notations) 

* CO- • <x> • -Offi- • ■„„ 

Full dot is the weak coupling vertex including all the diagrams of one sign, 
NN interaction block is the full block also of one sign diagrams. The lines 
are full Green functions with the widths. The most essential term is the one- 
loop diagram (see (45)), which is positive definite, and including the fcrmion 
width corresponds to the first term of the classical Langevin result, for details 
see [28]. Other diagrams represent interference terms due to rescatterings. In 
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sonic simplified representations (e.g., as we used within Fermi liquid theory) 
the 4-point functions (blocks of NN interaction of given sign diagrams) be- 
have like intermediate bosons (e.g. zero-sounds and dressed pions). In general 
it is not necessary to consider different quanta dealing instead with the full 
NN interaction (all diagrams of given sign). For particle propagation in an 
external field, e.g. infinitely heavy scattering centers (proper LPM effect), 
only the one-loop diagram remains, since one deals then with a genuine one- 
body problem. In the quasiclassical limit for fcrmions (with small occupations 
rip) all the diagrams given by first line of series (46) with arbitrary number 

of " h" lines are summed up leading to the diffusion result, for details 

see [28] . For small momenta q this leads to a suppression factor of the form 
C = uj 2 /(uj 2 + r 2 ), r x incorporates rescattering processes. In general case 
the total radiation rate is obtained by summation of all diagrams in (46). 
The value —ill h determines the gain term in the generalized kinetic equa- 
tion for G '"j see [67,68], that allows to use this method in non-equilibrium 
problems, like for description of neutrino transport in semi-transparent region 
of the neutrino-sphere of supernovae, as we may expect. 

In the QP limit diagrams 1, 2, 4 and 5 of (46) correspond to the MMU 
and MNB processes related to a single in-medium scattering of two fermionic 
QP and can be symbolically expressed as Feynman amplitude (47a) 

!©: 
(«) ■ "• : w *' 
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1, 1 (d) O— (47) 

" "*" * * 

The one-loop diagram in (46) is particular, since its QP approximant in many 

cases vanishes as we have mentioned. However the full one-loop includes 
QP graphs of the type (47b), which survive to the same order in r^ /u> c h 
as the other diagrams (therefore in [24] where QP picture was used this 
diagram was considered as allowed diagram). In QP series such a term is 
included in second diagram of (44) although beyond the QPA it is included 
as the proper self-energy insertion to the one-loop result, i.e. in first term 
(44) [28] . In fact it is positive definite and corresponds to the absolute square 
of the amplitude (47a). The other diagrams 2, 4 and 5 of (46) describe the 
interference of amplitude (47a) either with those amplitudes where the weak 
coupling quantum (Iv pair) couples to another leg or with one of the exchange 
diagrams. For neutral interactions diagram (46:2) is more important than 
diagram 4 while this behavior reverses for charge exchange interactions (the 
latter is important, e.g., for gluon radiation from quarks in QCD transport 
due to color exchange interactions) . Diagrams like 3 describe the interference 
terms due to further rescatterings of the source fermion with others as shown 
by (47c). Diagram (46:6) describes the production from intermediate states 
and relates to the Feynman graph (47d). For photons in the soft limit (to <C 
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Sf) this diagram (47d) gives a smaller contribution to the photon production 
rate than the diagram (47a), where the normal brcmsstrahlung contribution 
diverges like 1/u compared to the l/sp-vahie typical for the coupling to 
intermediate fermion lines. For vv brcmsstrahlung (47d) gives zero due to 
symmetry. However in some cases the process (47d) might be very important 
even in the soft limit. This is indeed the case for the MMU process considered 
above. Some of the diagrams which are not presented explicitly in (46) give 
more than two pieces, if being cut, so they never reduce to the Feynman 
amplitudes. However in the QPA they give zero contribution [28]. 

With Fp ~ tx 2 T 2 Jef for Fermi liquids, the criterion i> <C u> c h ~ T is 
satisfied for all thermal excitations As ~ T <C Ef/tt 2 . However with the 
application to soft radiation this concept is no longer justified. Indeed series 
of QP diagrams is not convergent in the soft limit and there is no hope 
to ever recover a reliable result by a finite number of QP diagrams for the 
production of soft quanta. With full Green functions, however, one obtains 
a description that uniformly covers both the soft (u> <C If) and the hard 
(uj 3> If) regimes. In the vicinity of g C7r the quantity TV being roughly 
estimated in [53,55] as TV oc 7T I Tm n /uj, and coherence effects come into 
play. 

In order to correct QP evaluations of different diagrams by the fermion 
width effects for soft radiating quanta one can simply multiply the QP results 
by different pre-factors [28]. E.g., comparing the one-loop result at non-zero 
rp with the first non-zero diagram in the QPA [Tp = in the fermion Green 
functions) we get 




= r„U')<-t f- ) , (48) 

QPA 
at small momentum q. For the next order diagrams we have 




GO- 




Ci(u){ -- + -4 O f- ,d(w)=w 



, , 2 r 2 
2 u> - l F 



QPA ^ + ^) 2 



/ \ + / N 




c\aM--\ D f ) . c oH = ^T7f • (49) 



where factors Co, C\, ... cure the defect of the QPA for soft uj. The factor Co 
complies with the replacement uj — > uj + i-Tp . A similar factor is observed in 
the diffusion result, where however the macroscopic relaxation rate r x enters, 
due to the resummation of all rescattering processes. 

Finally, we demonstrated how to calculate the rates of different reactions 
in dense equilibrium and non-equilibrium matter and compared the results 
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derived in closed diagram technique with those obtained in the standard 
technique of computing of the squared matrix elements. 

Concluding, the "nuclear medium cooling scenario" allows easily to 
achieve agreement with existing data. However there remains essential uncer- 
tainty in quantitative predictions due to a pure knowledge, especially, of the 
residual interaction treated above in an economical way within a phenomeno- 
logical Fermi liquid model which needs further essential improvements. 

Acknowledgement: Author appreciates the hospitality and support of 
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Abstract 

This review demonstrates that neutrino emission from the dense hadronic compo- 
nent in neutron stars is subject of strong modifications due to collective effects in 
the nuclear matter. With the most important in-medium processes incorporated in 
the cooling code an overall agreement with available soft X ray data can be eas- 
ily achieved. With these findings so called "standard" and "non-standard" cooling 
scenarios are replaced by one general "nuclear medium cooling scenario" which re- 
lates slow and rapid neutron star coolings to the star masses (interior densities). 
In-medium effects take important part also at early hot stage of neutron star evolu- 
tion decreasing the neutrino opacity for less massive and increasing for more massive 
neutron stars leading to a delay of the neutrino puis in the later case. An adequate 
formalism for calculation of neutrino radiation from strongly coupled matter is pre- 
sented that treats on equal footing one-nucleon and multiple-nucleon processes as 
well as resonance boson contributions. 



1 Introduction 



The EINSTEIN, EXOSAT and ROSAT observatories have measured surface temper- 
atures of certain neutron stars and put upper limits on the surface temperatures of 
some other neutron stars ( cf. [1-3] and further references therein). The data for 
some supernova remnants indicate rather slow cooling, while the data for several 
pulsars point to an essentially more rapid cooling. 

The physics of neutron star cooling is based on a number of ingredients, among which 
the neutrino emissivity of the high-density hadronic matter in the star's core plays 
a crucial role. Neutron star temperatures are such that, except first seconds-hours, 
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neutrinos/antineutrinos radiate energy directly from the star without subsequent 
collisions, since A„, \ 9 ^$> R, where X u , X p are the neutrino and antineutrino mean free 
paths and R is star's radius. In the so called "standard scenario" of the neutron star 
cooling (scenario for slow cooling) the most important channel up to temperatures 
T < (10 8 — 10 9 )K belongs to the modified Urea (MU) process nn — > npev. First 
estimates of its emissivity were done in refs [4,5]. Friman and Maxwell in refs. [6,7] 
recalculated the emissivity of this process in the model, where the nucleon-nucleon 
interaction was treated with the help of slightly modified free one-pion exchange 
(FOPE). Their important result for the emissivity, £„[FOPE], was proved to be by 
an order of magnitude larger than the previously obtained one. Namely the value 
£,,[FOPE] was used in various computer simulations resulting in the "standard sce- 
nario" of the cooling, e.g. cf. refs [8-10]. Besides the MU process, in the framework 
of the "standard scenario" numerical codes included also processes of the nucleon 
(neutron and proton) bremsstrahlung (NB) nn — >■ nnvv and np — >■ npvv, which 
lead to a smaller contribution to the emissivity then the MU, cf. refs [11,6]. Medium 
effects enter the above two-nucleon (MU and NB) rates mainly through the effec- 
tive mass of the nucleons which has a smooth density dependence. Therefore within 
FOPE model the density dependence of the rates is rather weak and the neutrino ra- 
diation from a neutron star depends only very weakly on its mass. This is the reason 
why the "standard scenario" based on the result of [6], though complying well with 
several slowly cooling pulsars, fails to explain the data of the more rapidly cooling 
ones. Also "standard scenario" included processes contributing to the emissivity in 
the neutron star crust which become to be important at a lower temperature. 

The non-standard scenario included so called exotics, associated with different types 
of direct Urea-like processes, i.e. the pion Urea (PU) [12] and kaon Urea (KU) [13,14] 
/3-decay processes and direct Urea (DU) on nucleons and hyperons [15] possible only 
in sufficiently dense interiors of rather massive neutron stars. The main difference in 
the cooling efficiency driven by the DU-like processes on one hand and the MU and 
NB processes on the other hand lies in the rather different phase spaces associated 
with these reactions. In the MU and NB case the available phase space is that of 
a two-fermion origin, while in the pion (kaon) /3-decay and DU on nucleons and 
hyperons it is that of a one-fermion origin. Critical density of pion condensation 
in neutron star matter is p C7T ~ (1 — 3)po depending on the type of condensation 
(neutral or charged) and the model, see [16-19]. Critical density of kaon condensation 
is p C K — (2 — 6)p depending on the type (K~ or K°, S or P wave) and the model, 
see [20,69]. Critical density for the DU process is p c u — (3 — 6)p depending on the 
model for the equation of state, see [15,18]. Most recent calculations [19] estimated 
critical density of neutral pion condensation as 2.5po an d for the charged one as 
1.7 po, whereas variational calculations [18] argued even for smaller critical densities 
(~ 1.3p for neutral pion condensation). On the other hand the equation of state of 
ref. [18] allows for DU process only at p > 5p . 

There is no bridge between "standard" and "non-standard" scenarios due to complete 
ignorance of in-medium modifications of NN interaction which allows for strong 
polarization of the soft modes (like virtual dressed pion and kaon modes serving 



a part of in-medium baryon-baryon interaction). Only due to increase of such a 
polarization with the baryon density mentioned condensates may appear and it 
seems thereby quite inconsistent to ignore these softening effects, i.e. precursors of 
the condensations for p < p C7r , and suddenly swith on the condensates for p > p C7T . 

Now let us basing on the results [22-25,16,26] briefly discuss a general "nuclear 
medium cooling scenario" which treats obvious caveats of two mentioned above 
scenarios. First of all one observes that in the nuclear matter many new reaction 
channels are opened up compared to the vacuum processes, see [22]. Generally speak- 
ing, standard Feynman technique fails to calculate in-medium reaction rates if the 
particle widths are important since there are no free particle asymptotic states in 
matter. Then summation of all perturbative Feynman diagrams where free Green 
functions are replaced by the in-medium ones leads to a double counting due to 
multiple repetitions of some processes (for an extensive discussion of this defect see 
[27]). This calls a formalism dealing with closed diagrams (integrated over all pos- 
sible in-medium particle states) with full non-equilibrium Green functions. Such a 
formalism was elaborated in [24] within quasiparticle approximation for nucleons 
and was called "optic theorem formalism (OTF) in non- equilibrium diagram tech- 
nique". It was demonstrated that both standard calculation of the rates via squared 
reaction matrix elements and calculation using OTF coincide within quasiparticle 
approximation picture for the fermions. In [28] the formalism was generalized to in- 
clude arbitrary particle widths effects. The latter formalism treats on equal footing 
one-nucleon and multiple-nucleon processes as well as resonance reaction contribu- 
tions of the boson origin, as processes with participation of zero-sounds and reactions 
on the boson condensates. Each diagram in the series with full Green functions is 
free from the infrared divergences. Both, the correct quasiparticle and quasiclassical 
limits are recovered. 

Except for very early stage of neutron star evolution (seconds - hours) typical av- 
eraged lepton energy (> T) is larger then the nucleon particle width T N ~ T 2 /e FN 
and the nucleons can be treated within the quasiparticle approximation. This ob- 
servation much simplifies the consideration since one can use an intuitive way of 
separation of the processes according to the available phase space. The one-nucleon 
processes have the largest emissivity (if they are not forbidden by energy-momentum 
conservations), then two-nucleon processes come into play, etc. 

In the temperature interval T c < T < T opac (T c is typical temperature for the nu- 
cleon pairing and T opac is typical temperature at which neutrino/antineutrino mean 
free path \ v j\ v becomes to be of the order of the star radius R) the neutrino 
emission is dominated by the medium modified Urea (MMU) and medium nucleon 
bremsstrahlung (MNB) processes since one-nucleon processes like DU, PU and KU, 
are forbidden up to sufficiently large nucleon density. Corresponding diagrams of 
MMU and MNB are schematically shown in Fig. 1. Refs [22-25,16] considered NN 
interaction within Fermi liquid Landau-Migdal approach. They incorporated the 
softening of the medium one-pion exchange (MOPE) mode, other medium polariza- 
tion effects, like nucleon-nucleon correlations in the vertices, renormalization of the 
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Fig. 1. Antineutrino emission from a nucleoli leg (left graph) and from intermediate scat- 
tering states (right). Full dot includes weak coupling vertex renormalization. 
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Fig. 2. Neutrino-antineutrino emission from Cooper pair-breaking (left graph) and 
pair- formation processes (right graph). 

local part of NN interaction by the loops, as well as the possibility of the neutrino 
emission from the intermediate reaction states and resonance DU-like reactions going 
on zero sounds and the boson condensates. Refs [23,25,16] have demonstrated that 
for p > pq second diagram of Fig. 1 gives the main contribution to the emissivity of 
MMU process rather than the first one which contribution has been evaluated much 
earlier in the framework of FOPE model in ref. [6]. This fact essentially modifies 
the absolute value as well as the density dependence of the nn — >■ npev process rate 
which becomes to be very sharp. Thereby, for stars of masses larger than the sollar 
mass the resulting emissivities were proved to be substantially larger than those val- 
ues calculated in FOPE model. With increase of the star mass (central density) pion 
mode continues to soften and MMU and MNB rates still increase. At p > p C7r pion 
condensation begins to contribute. Actually, the condensate droplets exist already 
at a smaller density in the mixed phase appearing in systems having more than 
one conserved charge [29,30]. At T > T me i t , where T me i t is the melting temperature, 
roughly ~ several MeV, the mixed phase is in liquid state and PU processes on 
independent condensate droplets are possible. At T < T me i t condensate droplets are 
placed in a crystalline lattice that substantially suppresses corresponding neutrino 
processes. 

Ref. [12] considered the reaction channel n — > pit'ev, whereas refs [22,23] included 
other possible pion 7r + , tt ± , 7r° condensate processes with charged and neutral cur- 
rents (e.g., like nn^ — >■ puu and nn® — >■ nvv) as well as analogous resonance reactions 
going on zero sounds which are also possible at p < p C7r . Due to AW correlations all 
pion condensate rates are significantly suppressed (by factors ~ 10 - 100 compared 
to first estimate [12], see [31,22,23]. At p ~ p C7T both MMU and PU processes are of 
the same order of magnitude showing a smooth transition to higher densities (star 
masses) in the "nuclear medium cooling scenario" being absent in the "standard" 
and "non-standard" scenarios. 

For T < T c the reactions of neutrino pair radiation from superfluid nucleon pair 
breaking and formation (NPBF) shown in Fig. 2 become to be dominant processes. 
The neutron pair breaking and formation (nPBF) process for the case of the neutron 



liSo pairing was first calculated in ref. [32] using standard Bogolyubov technique. 
Later this process was independently calculated in [24,25] as demonstration of ef- 
ficiency of OTF within the closed non-equilibrium diagram technique developed 
there. Moreover ref. [24] calculated emissivity of the corresponding process on pro- 
ton (pPFB) taking into account strong coupling pvv vertex renormalization (see 
first diagram eq. (20) below). It resulted in a hundred times enhancement in the 
medium of pPFB emissivity compared to that estimated with the vacuum vertex, 
leading to that both nPFB and pPFB emissivities were proved to be of the same or- 
der of magnitude. Emissivities of NPBF processes have the same suppression factor 
~ exp(-2A/T) as MU, NB, MMU and MNB processes at T < T c but compared to 
the latter the NPBF processes have a large one nucleon phase space volume. Ref. 
[35] included peculiarities of 3P2 pairing. In the 3P2G m J |= 2) case, where mj is 
projection of the total pair momentum onto quantization axis, exponential suppres- 
sion of specific heat and the emissivity is replaced by only a power law suppression 
since the gap vanishes at the Fermi sphere poles, as was first claimed in ref. [24] and 
then included in refs [36,33] which overlooked this more earlier suggestion. 

At this instance I would like to do a historic remark relating to estimation of the 
NPBF processes since in some refs (see [34,33]) was expressed a surprize why these 
processes were not on a market during many years. First work [32], although found 
correct analytic expression for 1S pairing of neutrons, numerically underestimated 
the emissivity by more than an order of magnitude (by factors 20 and 50 for S 
and P paring cases respectively). Also there was no statement on the dominance 
of this process in the cooling scenario (over the MU). Asymptotic behaviour of the 
emissivity e[nPBF] ~ 10 20 T 9 7 exp(-2A/T) for T < A, T 9 = T/10 9 K, as follows 
from their expression (lb) of [32] and from their rough asymptotic estimate of the 
integral (see below (13b)), shows up nor a large one nucleon phase space factor 
(~ 10 28 ) nor appropriate temperature behaviour. Namely this underestimation of 
the rate in [32] (we again point out that analytic expression (la) is quite correct) 
and absence of mentioning of a possible dominance of the process over MU, became 
a reason that this important result was overlooked during many years. Ref. [24] 
overlooked sign of anomalous diagram for S pairing and included ~ g\ term, where 
qa is axial- vector coupling constant, which should absent in the case of S pairing 
but being the most important contribution for the 3P2 pairing case. In spite of 
that, a reasonable numerical estimate was presented valid for both S and P pairing 
including NN correlation effects into consideration. Uncertainty of this numerical 
value can be estimated by a factor (0.5-1) which is in any case allowed by variation of 
not too well known correlation factors. Namely this estimate was then incorporated 
within the cooling code in [26]. Correct asymptotic behaviour of the emissivity is 
e[nPFB,pPFB] ~ 10 28 (A/MeV) 7 (T/A) 1 / 2 exp(-2A/T) for T < A, that showed 
up a huge one-nucleon phase space factor and very moderate T dependence of the 
pre-factor. Thereby the value of the rate was related in [24] to the value of the pairing 
gap. The possibility of the dominant role of this process (even compared to enhanced 
MMU rather than only MU) was unambiguously stressed. Unfortunately ref. [24] had 
a number of obvious misprints which were partially corrected in subsequent papers 
[25,16,35]. Although it does not excuse the authors these misprints can be easily 



treated by an attentive reader. Ref. [26] was the first that found and quated the 
previous result [32]. 

Ref. [26] incorporated most important in-medium effects in the cooling code, among 
them nPBF and pPBF as equally important processes. However importance of pPBF 
process was then overlooked in the subsequent papers where emissivity of this pro- 
cess was several times incorrectly reproduced. Subsequent works [34,35] supported 
conclusion of [26] on importance of NPBF processes as governing the cooling scenario 
at T < T c . 

The medium modifications of all the above mentioned rates result in a pronounced 
density dependence (for NPBF processes mainly via dependence of the pairing gaps 
on the density and the NN correlation factors), which links the cooling behavior of 
a neutron star decisively to its mass. As the result, the above mentioned medium 
modifications lead to a more rapid cooling than obtained in the "standard scenario". 
Hence they provide a possible explanation for the observed deviations of some of 
the pulsar temperatures from the "standard" cooling. Particularly, they provide a 
smooth transition from "standard" to "non-standard" cooling for increasing cen- 
tral star densities, i.e., star masses. Thus by means of taking into account of most 
important in-medium effects in the reaction rates one is indeed able to achieve an 
appropriate agreement with both the high as well as the low observed pulsar temper- 
atures that leads to the new "nuclear medium cooling scenario". Using a collection 
of modern equations of state for nuclear matter in its ground state, which covered 
both relativistic as well as non-relativistic models, ref. [26] also has demonstrated a 
relative robustness of the new cooling mechanisms against variations in the equation 
of state of dense neutron star matter. 

At initial stage (T > T opac ) the newly formed hot neutron star is opaque for neutri- 
nos/antineutrinos. Within FOPE model the value T opac was estimated in [6]. Elas- 
tic scatterings were included in ref. [37,38] and pion condensation effect on the 
opacity was discussed in ref. [39]. Medium effects drammatically affect the neu- 
trino/antineutrino mean free paths, \ v („) ex 1/ | M | 2 , where M is the reaction 
matrix element. Thereby, one-nucleon elastic scattering processes, like Nv — > Nv, 
are suppressed by AW correlations [24,16,40,41]. However neutrino/antineutrino ab- 
sorption in two-nucleon MMU and MNB processes increases with the density (since 
| M | 2 of MMU and MNB processes increase with the density) [23,24,16]. Thus 
more massive neutron stars are opaque for neutrinos/antineutrinos up to lower tem- 
peratures that also results in a delay of neutrino/antineutrino puis. Within the 
quasiparticle approximation for the nucleons the value T opac was estimated with 
taking into account medium effects in refs [23,16]. Refs [42,43,16] considered pos- 
sible consequences of such a delay for supernova explosions. On the other hand 
at T > (1 — 2)MeV one should take care of the neutrino/antineutrino radiation 
in multiple AW scatterings (Landau-Pomeranchuk-Migdal effect) when averaged 
neutrino-antineutrino energy, uj v9 ~ several T, becomes to be smaller than the nu- 
cleon width T N [28]. Numerical evaluations of T N in application to MNB processes 
were done in [76] and [77] within the Bruckner scheme and using Bethe-Salpeter 



equation for nucleons. The Landau-Pomeranchuk-Migdal effect suppresses the rates 
of the neutrino elastic scattering processes on nucleons and also it suppresses MNB 
rates. For neutron stars of rather low mass ( ^ -^q) the suppression of the rates of 
neutral current processes due to NN correlations and the multiple collision coher- 
ence effect prevail the enhancement due to the pion softening, and for sufficiently 
massive neutron stars (> 1.4Mq) the enhancement prevails the suppression. MMU 
emissivity remains to be almost unaffected by the Landau-Pomeranchuk-Migdal 
effect since averaged ve energy ~ p Fe is rather large ^> T N . 

The paper is organized as follows. Sect. 2 dicusses basic ideas of the Landau-Migdal 
Fermi liquid approach. The AW interaction amplitude is constructed with an explicit 
treatment of long-ranged soft pion degree of freedom and vertex renormalizations 
due to AW correlations. The meaning of the pion softenning effect is clarified and a 
comparison for AW cross sections evaluated in MOPE and FOPE models is given. 
Then the renormalization of weak interaction in neutron star matter is performed. 
Sect. 3 discusses cooling of non-superfluid neutron stars (T opac > T > T c ). We start 
with a demonstration that it is theoretically inconsistent to describe AW interac- 
tion with the help of vacuum cross section at the same time considering pions as 
freely propagating particles in matter. This inconsistency will lead us to a problem 
which will be then resolved within a consistent Fermi liquid approach. The com- 
parison of emissivities of MMU and MU processes shows a significant enhancement 
of in-medium rates. Then we discuss DU-like processes and demonstrate a suppres- 
sion effect due to vertex renormalizations. Demonstration of the role of in-medium 
mechanisms in the calculation of the cooling evolution of neutron stars within a real- 
istic cooling code is given in next Section. Then we consider influence of in-medium 
effects on the neutrino mean free path at initial stage of neutron star cooling. Essen- 
tial role of multiple AW collisions (Landau-Pomeranchuk-Migdal effect) is shown. 
Sect. 4 presents OTF in non-equilibrium closed diagram technique in the framework 
of quasiparticle approximation for the nucleons and then beyond the quasiparticle 
approximation incorporating particle width effects. 



2 Fermi Liquid Behaviour of Baryon Matter 

2.1 A general strategy - separation of soft and hard modes 



We will restrict ourselves by consideration of the nuclear matter at densities less 
than about 5po, where po — 0.17 fm~ 3 is the normal nuclear matter density, and 
at typical excitation energies (temperatures) which are very low compared to those 
essentially involving quark-gluon degrees of freedom (~ l/r\, where r^ ~ (0.2 — 0.3) 
fm is the confinement radius). Then, we may deal with baryon-meson degrees of 
freedom rather than with the quark-gluon ones. Thus, we suppose that baryons 
interact via meson (pion, sigma, rho, omega, etc.) exchanges. 



Different processes involved in nuclear phenomena occuring on hadron level are 
characterized by the corresponding space-time and energy-momentum scales. Of the 
order of the confinement radius r^ are scales associated with heavy mesons such as a 
and Q mesons (r ~ r a ~ l/m„, r ~ r^ ~ 1/mo) and with nucleons, r^ ~ 1/mjv- We 
will treat these scales as short-ranged. Typical low momenta in our consideration 
will be k ~ m^ and k ~ p F ~ l/(2m 7r ), where p^ is Fermi momentum of nucleon, 
that corresponds to the length scales r ~ 1/m^ and r ~ 1/pf, respectively. These 
latter scales will be treated as rather long-ranged microscopic scales. Time scales, 
which correspond to low energies are ~ 1/ep (ep is the Fermi energy of nucleons, 
~ l/m v and ~ 1/uja ( w a = rnA — t^n — 2.1m v is the energy gap for the excitation 
of A isobars). A isobars should therefore be treated as long-ranged excitations in 
our scales. Besides they have a large spin-isospin degeneracy factor (16) compared 
to that for nucleons (4) and more strongly coupled to the pions (f V NA — 2.151 /m n 
whereas f^N — V 77 ^)- Especial role here belongs to the pion mode. The saturation 
nuclear density is p ~ 0.48m^, i.e. of order of unit in these natural interaction units. 
Already from this simple observation one can conclude first that pion takes essential 
part in gluing of nuclear matter and second, that the nucleon-nucleon interaction 
and, in particular, long-range part of this interaction at p > p should essentially 
differ from those in vacuum. Thereby, there is, actually, no hope that perturbative 
schemes could help in construction of a proper description of nucleon matter for 
p > po- In such a situation, one should expect strong in-medium effects, especially 
for pions, being the lightest among hadrons (m v /m N ~ 1/7). Minimal energy scale 
related to this long-range part of NN interaction is the effective pion gap lo < m n . 
The value u 2 replaces the corresponding value ~ m^ + p\ characterizing FOPE 
contribution to NN interaction, cf. [16]. The energy scale uj arises due to a peculier 
collective pion softening effect appearing for nucleon densities p > p c \ ~ (0.5 — 0.7) po 
due to a strong P wave ttN attraction discovered by Migdal [44,67,16]. Presence of 
this soft mode becomes of prime importance. As ions in solids interact mainly by 
phonons, i.e. by soft collective modes being absent in vacuum, as nucleons in a dense 
nucleon matter interact via the collective modes the dominant among which is the 
softest one. The latter is the particle-hole mode with the pion quantum numbers. 



2.2 Description of NN interaction 



At temperatures of our interest (<C ep n ) neutrons are only slightly excited above 
their Fermi sea and all the processes occur in a narrow vicinity of tp n . In such a 
situation Fermi liquid approach seems to be the most efficient one. Within this ap- 
proach the long-ranged diagrams are treated explicitly whereas short-scale diagrams 
are treated as local quantities given by phenomenological so called Landau-Migdal 
parameters. Thus using argumentation of Fermi liquid theory [46,47,67,16] the re- 
tarded A^A^ interaction amplitude is presented as follows 



X 





(1) 



where 




X + X 



(2) 





The solid line presents the nucleon, whereas double-line, the A isobar. The empty- 
dotted line corresponds to the exchange of the free pion with inclusion of the con- 
tributions of the residual S wave ttNN interaction and 7T7t scattering, i.e. the ir- 
reducible residual interaction to the nucleon particle-holes and delta-nucleon holes. 
The full particle-hole, delta-nucleon hole and pion irreducible block in eq. (1) is by its 
construction essentially more local then contributions given by explicitly presented 
graphs. Thereby, it is parameterized with the help of the so called Landau-Migdal 
parameters. In the standard Landau Fermi liquid theory fermions are supposed 
to be at their Fermi surface and the Landau parameters are expanded in Legen- 
dre polinomials in the angle between fermionic momenta. Luckily, only zero and 
first harmounics enter physical quantities. The momentum dependence of nuclear 
fources is not as pronounced and one has no need to do this expansion. Then these 
parameters, i.e. f nn , f np and g nn , g np in scalar and spin channels respectively, can 
be considered as slightly momentum dependent quantities. They can be evaluated 
as functions of density, neutron and proton concentrations, energy and momentum. 

The part of interaction involving A is analogously constructed 

(3) 



The main part of the NA interaction is due to the pion exchange. Although informa- 
tion on local part of the iVA interaction is rather scarce one can conclude [16,19] that 
the corresponding Landau-Migdal parameters are essentially smaller then those for 
NN interaction. Besides, at small transferred energies u <C m, the A-nucleon hole 
contribution is a smooth function of u and k in difference with the nucleon-nucleon 
hole contribution. Therefore and for the sake of simplicity we neglect the first graph 
in r.h.s. of eq. (3). 

Straightforward resummation of (1) in neutral channel yields 

a (5 

r^=y^= C (^ + Z%a x ■ a 2 ) + f vN T« (ffi ■ k) (a 2 • k) , (4) 

a (5 

Faf3 = fapl(faf)), %nn = g n nl{9nn) , %np = 9npl(9nn) , <*, P = (", p) , 
Tnn = l 2 (9nn)D* , T* p = -J pp j(g nn )D^ , T pp = J^pD^o, 

1 -\x) = l-2C Q xA* n , %p = (l-4C gA^ n ) 1 (g nn ), (5) 

fnn = fpp = /+/', fnp = f-f, 9nn = 9 P p = 9+9', and g np = g-g', C is dimensional 
normalization factor which is usually taken to be Co = 7r 2 /[m N p F (p )] — 300 MeV- 



fm 3 ~ 0.77m v 2 , D^ is the full retarded Green function of 7r°, A a p is the correspond- 
ing (nn) nucleon-nucleon hole loop 

A«p = (^) (6) 



a- 1 



and we for simplicity neglect proton hole contributions due to a small concentration 
of protons in the neutron star matter. 



Resummation of (1) in the charged channel yields 
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f * = V = C (fl + Z%oi ■ B 2 ) + fl N f^ ■ k)(a 2 ■ k) 
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•^ = 2/'7(A Z«=2gft(gf) t f£ = ?tf)D?- , T\x) = 1 - 4C xA^ p . (8) 

The Landau-Migdal parameters are rather unknown for isotopically asymmetric nu- 
clear matter and for p> p . Althogh some evaluations of these quantities have been 
done, much work is still needed to get covincing results. Therefore for estimates we 
will use the values extracted from atomic nucleus experiments. Using argumenta- 
tion of a relative locality of these quantities we will suppose the Landau-Migdal 
parameters to be independent on the density for p > po- One then can expect that 
the most uncertain will be the value of the scalar constant / due to essential role of 
the medium-heavy a meson in this channel. But this parameter does not enter the 
tensor fource channel being most important in our case. Unfortunately, there are 
also essential uncertainties in numerical values of some of the Landau-Migdal pa- 
rameters even for atomic nuclei. These uncertainties are, mainly, due to attempts to 
get the best fit to experimental data in each concrete case slightly modifying param- 
eterization used for the residual part of the A^A^ interaction. E.g., calculations [47] 
gave / ~ 0.25, /' ~ 1, g ~ 0.5, g' ~ 1 whereas refs [48-50], including quasiparticle 
renormalization pre-factors, derived the values / ~ 0, /' ~ 0.5 — 0.6, g ~ 0.05 ± 0.1, 
0'~ 1.1 ±0.1. 

Typical energies and momenta entering A^A^ interaction of our interest are lo ~ and 
k ~ p Fn . Then rough estimation yields 7~ 1 (g / , u) — 0, k ~ p Fn , p = p ) ~ 0.3 — 0.45. 
For u) = k ~ T typical for the weak processes with participation of vv one has 
7" 1 (y,a> ~ k~T,p = po) ~ 0.7 -0.9. 



2. 3 Virtual Pion Mode of NN Interaction. 

Resummation of the diagrams yields the following Dyson equation for pions 

= ~ +~o +0 + ^m^ 
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(9) 



The ttNA full-dot-vertex includes a background correction due to the presence of 
the higher resonances, U r r f s is the residual retarded pion self-energy that includes 
the contribution of all the diagrams which are not presented explicitly in eq. (9). 
The full vertex takes into account NN correlations 
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Due to that the nucleon particle-hole part of II^o is oc j(g nn ) and the nucleon 
particle-hole part of U n ± is oc 'y(g'). The value of the A^A" interaction in the pion 
channel is determined by the full pion propagator at small lo and k ~ pp n , i.e. 
by the quantity L0 2 (k) = — (D^)~ 1 (lo = 0,k,/j, v ). Typical momenta of our interest 
are k ~ pp n . Indeed the momenta entering the A^A" interaction in MU and MMU 
processes are k = pp n , the momenta governing the MNB are k = k [23] where the 
value k = k ~ (0.9 — l)pF n corresponds to the minimum of u) 2 (k). The quantity 
lo = Lo(ko) has the meaning of the effective pion gap. It is different for n° and 
for 71"^ since neutral and charged channels are characterized by different diagrams 
permitted by charge conservation, thus also depending on the value of the pion 
chemical potential, //+ ^ jjl~ ^ 0, //° = 0. For T <C ep n , £f p these chemical potentials 
are p~ = (i e = tp n — £f p , [^ = e Fn, as follows from equilibrium conditions for the 
reactions n — > pn~ , p — >■ nn + and n — > pev. 

Change the sign of lo 2 symbolizes the pion condensate phase transition. Typical 
density behaviour of lo 2 is shown in Fig. 3. At low densities (p < (0.5 — 0.7)po), 
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Fig. 3. Effective pion gap (for n^o 
RelT? is small and one has lo 2 



mz 



0) versus baryon density, see ref. [16]. 

/j, 2 . For such densities the value u 2 (p Fn ) 



essentially deviates from m\ — //| tending to the value m\ 
density limit. Approximately one has fi n - = pp e = pp p = PFn/Z. 



Pin 



fi^ only in small 



The dashed line describes the case where 7T7t fluctuations are artificially switched off 
and the phase transition turns out to be of the second order. At the critical point of 
the pion condensation (p = p C7T ) this value of lo 2 with switched off mr fluctuations 
changes its sign. In reality the mr fluctuations are significant in the vicinity of the 
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critical point and there occurs the first-order phase transition to the inhomogeneous 
pion-condensate state [51,59,53]. Thereby there are two branches (solid curves in 
Fig. 3) with positive and respectively negative value for uj 2 . Calculations of ref. [53] 
demonstrated that at p > p C7T the free energy of the state with u 2 > 0, where the pion 
mean field is zero, becomes larger than that of the corresponding state with ui 2 < 
and a finite mean field. Therefore at p > p c the state with u 2 > is metastable and 
the state with u 2 < and the pion mean field ip n ^ becomes the ground state. 

The quantity ui 2 demonstrates how much the virtual (particle-hole) mode with pion 
quantum numbers is softened at given density. The ratio a = D 7r [med.]/.D 7r [vac.] ~ 6 
for p = po and for isospin symmetric nuclear matter. However this essential so called 
"pion softening" does not significantly enhance the NN scattering cross section due 
to a simultaneous suppression of ttNN vertex by NN correlations. Estimation of 
the rate of the NN cross sections calculated via FOPE and MOPE yields 

_ q[FOPE] i*(g',u}~0,k~p FN )(ml+p 2 FN ) 

a[MOPE] " £ 4 ' [ ' 



For p = po one has jig'^co ~ 0, k ~ Pfn) — 0.3 — 0.4, lu 2 ~ 0.8, p F N — 2m,,- and 
R& 1, whereas for p = 2p with our rough estimates one already gets R ~ 10. 

As follows from numerical estimates of different 7 factors entering (4) and (7) the 
most important contribution to A^A^ interaction for p > p is given by MOPE 



X 



(12) 



whether this channel (T oc (e?i • k)(a 2 ■ k)) of the reaction is not forbidden or 
suppressed by some specific reasons like symmetry, small momentum transfer, etc. 
Thus instead of FOPE, as the model of AW interaction which has been used in 
[6] in calculation of the emissivities of two-nucleon reactions within the "standard 
scenario" of neutron star cooling one needs to use the full AW interaction given by 
eqs (4) and (7) or, simplifying, approximated by its MOPE part. 

Some doubts were expressed in connection to the role of the p meson part of tensor 
fources. Ref. [6] estimated p meson contribution to MU emissivity arguing for a 
suppression of FOPE contribution by factor of 2. Let us estimate how much is 
essential a contribution of the p meson to the tensor fources. Using identity (e?i x 
k){&2 x k) = k 2 ai&2 — ipik) ■ (a 2 k) the p exchange interaction can be cast in the 
form 



r r n = (_fp\ I k (gifW (a 1 fc)(q 2 A:)7 2 

NlN2 \mj \[uj 2 -m 2 p -k 2 -Uf} [co 2 -m 2 p -k 2 -Uf] r 
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We may omit contribution Rell^ in the Green function since Rell^ <$; m 2 for p ~ m\ 
under consideration; 7 factor is the same as for the corresponding pion (charged 
or neutral). First term is supposed to be included in phenomenological value of 
the corresponding Landau-Migdal parameter leading to its momentum dependence. 
Due to that the value g' gets a 30% decrease (rather then an increase discussed in 
some works) with the momentum at k = p F N compared to the corresponding value 
at k = 0, see [56]. Thus one can think that g'(0) > g'(pFN)- Second term can be 
droped as small ( ^ 1/30 of MOPE contribution at p = p ). 

Another warry [57,58] was expressed in connection with experimental quasielastic 
polarization transfer experiment at LAMPF, EMC experiment and the Drell Yan 
experiment at Fermi Lab which did not observe an expected in several works a 
pronounced pion excess in nuclei. First optimistic estimates demonstrated the pion 
excess to be at (15-30)% level that was ruled out by different model analysis of men- 
tioned experimental results which included artificially enhanced pion contribution 
and suppressed other important contributions. There were also many rather naive 
speculations which expected the excess to essentially grow in the vicinity of the pion 
condensation point for (T = 0). 2 

First let us to briefly discuss the physical origin of the effect. Let us estimate a 
contribution of the virtual pion sea (region of small energies uj < kvp) to the full 
sum rule. It is given by 



S(kv F ) = 


kv F 

r , duj 2 

= J ^ A 'T, = Tp 




S(oo) = 


= 1, l3= m » k F N , 



u 2 (k) ((5kv F \ 

kv F — arctan ^ 9/; , 

\u 2 {k)J 



<1, 



(14) 



m* N is the nonrelativistic effective mass given by the relation {de(p)/dp)\ PF = pp/m* . 
The quantity S(kvp) having no singularity even at LO 2 (k ) — >■ is, thus, rather in- 
sensitive to the pion softening effect. Therefore it seems clear that the mentioned 
experimental analysis is not critical to the value ui 2 (ko). This point was missed in 
previous analysis of the pion excess. The physical reason of a non-zero contribu- 
tion of virtual pions to the sum-rule at u < kv F is clear. It is just the well known 
Landau damping associated with the possibility of the virtual pion decay to the 
nucleon particle and the hole, i.e. just with the Pauli blocking. As known from ttN 
scattering experiments, pions interact with nucleons. Thus must be a contribution 
to the sum-rule. There is no other way out. Therefore on a question "Where virtual 
pions are?" we would still suggest an old fashion answer: "They are hidden inside 
the matter". The question which remains is only "What is the actual value of the 
pion excess?" To properly answer it one needs besides the pionic contribution to 
ellaborate all relevant specific contributions related to the phenomenon analysed 
in the given concrete experiment. E.g., quasielastic polarization transfer data are 



For T/0 this excess indeed essentially grows in the vicinity of the critical point 
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reasonably fitted with a slightly increased value of the Landau-Migdal parameter 
g' and with inclusion of S wave repulsion at small energies [16,60,61], EMC exper- 
iment is reproduced with inclusion of the soft pion mode and other relevant effects 
[62], Drell-Yan production is described with taking into account of nucleon binding 
[63] and initial state interactions of fast partons [64] with the pion excess included 
(although in [64] the later is taken in a rather naive model where the pion excess is 
artificially increased) . 



2.4 Renormalization of the Weak Interaction. 



The full weak coupling vertex that takes into account NN correlations is determined 
by eq. (15) where now the wavy line should be replaced by the lepton pair. 



><=><+ ^< 



(15) 



Thus for the vertex of our interest, Ni — >■ iV 2 ^, we obtain 



V p = -£= [t(/% - 9aWM , (16) 



V2 
for the (5 decay and 



G r ,. ,, , ,7^1 „ G 



V nn = 7—7= \l(fnn)k ~ 9 Al{.9nn)la\ , V pp = —j= [n pp l - g A lpplo\ , (17) 



k pp = c v - 2f np C A nn 'j(f nn ), jpp = (1 - 4gC A nn ) j{g nn ), (18) 

for the processes on the neutral currents N\N 2 — >■ N\N2Vv. Here G ~ 1.17 • 
1CT 5 GeV~ 2 is the Fermi weak interaction constant, cy = 1 — 4sin 2 9w, sin 2 9\y — 
0.23, qa — 1-26 is the axial- vector interaction coupling constant, and l^ = ^(#1)7^(1 — 
75)^(^2) is the lepton current. The tensor fource contribution ~ q 2 is small for typ- 
ical I q |~ T or pp e , and for the sake of simplicity is ommited. 

The 7 factors renormalize the corresponding vacuum vertices. These factors are 
essentially different for different processes involved. The matrix elements of the neu- 
trino/antineutrino scattering processes Nv — > Nv and of MNB, NN — > NNvv, 
behave differently in dependence whether N = n or N = p m the weak coupling 
vertex. Vertices 



(19) 



V 


V 




V 




1 


N -^ 


1$° 


V 


N 
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N 



are modified by the correlation factors (5) and (8). For N = n these are ^j{g nni uj = 
q ~ T) and 7(/ rm ,tJ = q ~ T) which are ~ 0.7 — 0.9 leading to the suppression of 
the corresponding cross sections by 0.5 — 0.8. Renormalization of the proton vertex 
is determined by the processes 





(20) 



being forbidden in vacuum. For the systems with lSo proton-proton pairing, oc g\ 
contribution to the squared matrix element from the first diagram is compensated 
by the corresponding contribution of the diagram with anomalous Green functions 
of protons. The vector current term is oc Cy in vacuum. Thereby the corresponding 
vertices with participation of proton (N = p) in medium are enhanced compared 
to the small vacuum value (ex c v ~ 0.06) leading to the enhancement of the cross 
sections by ~ 10 2 times. Also contribution of the same order of magnitude may come 
from the intermediate in-medium photon j m whose propagator contains 1/w 2 / ~ 
1/e 2 (or 1/m 2 ~ 1/e 2 in superconducting proton medium) that compensates small 
e 2 factor from electromagnetic vertices [54]. Other permitted intermediate states 
like processes with pp~ l and with the pion are suppressed, the former by a small 
proton density factor and the latter by small q 2 ~ T 2 pre-factors. First diagram was 
considered already in [24] where the pPBF process was suggested, and then in [26], 
and it was shown that nPBF and nPBF processes may give contributions of the 
same order of magnitude. Several subsequent papers overlooked these results and 
rediscovered the pPBF process considering its rate as small since they ignored the 
nucleon and the electron correlation effects. Specific contribution of second diagram 
was recently discussed in [55] in association with pPBF process. First diagram was 
overlooked. The sub-process 



v\ 




(21) 



K- 



of v scattering on nucleus with excitation of in-medium K~ is drammatically mod- 
ified by in-medium effects [27]. The /i + production cross section is larger than the 
positron one by a factor fa 4 in contrast to expectation based upon the vacuum 
branching rations of a kaon decay T(K~ — >■ e~ + h , e )/T(K~ — > jjT + z^) ~ 2.5 x 10~ 5 . 
For the bare weak interaction the squared matrix element of the reaction T>i — >■ 
K~ + l + would be \M.\ 2 ex m 2 [m 2 — (k ■ k)]/2, that explains a strong enhance- 
ment of muon processes in vacuum compared to the positron ones. In medium the 
weak kaon current is dramatically modified due to the mixture of kaons with the 
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A-particle-proton hole states carrying the same quantum numbers 



~4 = «~f + — *CX w 



As the result, the squared matrix element of the reaction does not possess strong 
dependence on the lepton mass. Therefore, in medium squared matrix elements for 
positrons and muons turn out to be of the same order of magnitude. 

Thus we see that in dependence of what reaction channel is considered in-medium 
effects may or strongly enhance the reaction rate under consideration or substantially 
suppress it. 



2. 5 Inconsistency of FOPE Model for the Description of NN Interaction in Dense 

Matter 



Since FOPE model became the base of the "standard scenario" for cooling simula- 
tions we would like first to demonstrate principal inconsistencies of the model for 
the description of interactions in dense (p > po) baryon medium. The only diagram 
existing in FOPE model which contributes to the neutrino emissivity of two-nucleon 
processes is as follows 




(23) 



Dots symbolize FOPE. This is first available Born approximation diagram, i.e. sec- 
ond order perturbative contribution in f 2 N . Then in order to be theoretically consis- 
tent one should use perturbation theory up to the very same second order in f nN for 
all the quantities. E.g., pion spectrum is determined by pion polarization operator 

2 



expanded up to the very same second order in f 2 N 

u 2 = ml + k 2 + n Q {u,k,p), U°(u,k,p)= /7 T^V (24) 

The value U°(u,k,p) is easily calculated containing no any uncertain parameters. 
For the values u — >■ and k ~ pp of our interest one gets 



2m N ppk 2 f 2 N ^ n a m 2 N kflj 

r2 



n u ~ -« + *A)W, «o = \ nN > o, A, = -^^- > o. (25) 

7T Z TV 



Replacing this value to eq. (24) we obtain a solution with iu < already for p > 
0.2po that would mean appearance of the pion condensation. Indeed, the mean field 
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begins to increase with time passage (p ~ exp(— iut) ~ exp(ctt//5) untill repulsive 
7T7T interaction will not stop its growth. But it is experimentally proven that there 
is no pion condensation in atomic nuclei, i.e. even at p = p . The puzzle is solved as 
follows. FOPE model does not work for such densities. One should replace FOPE by 
the full NN interaction given by (4), (7). Essential part of this interaction is given 
by MOPE with vertices corrected by NN correlations. The nucleon-nucleon hole 
part of the pion polarization operator is also corrected by NN correlations. Thus 



o- 



iT>,M7(<?>,M (26) 



being suppressed by the factor j(g', uj = 0, k ~ p F , p ~ p ) ~ 0.3 — 0.4. Final solution 
of the dispersion relation (24) now with full n instead of n yields iuj > for p = p 
whereas the solution with iu < 0, which demonstrates the possibility of the pion 
condensation, appears only for p > p C7r > p . 



3 Neutrino Cooling of Neutron Stars 

3.1 Emissivity of MMU process. Comparison between MMU and MU 



Going from smaller to larger densities we have the following picture (for T opac > 
T > T c ). Since DU-like processes are forbidden up to sufficiently high density, the 
main contribution for T > T c comes from MMU processes schematically given by 
two diagrams of Fig.l. Contribution of MNB reactions is somewhat smaller except 
a narrow vicinity of the tt condensate point [23]. 

For densities p < 0.5 — 0.7 p the main part of the A^A" interaction amplitude is given 
by the residual A^ A" interaction that was approximated above by Landau-Migdal pa- 
rameters rather than by MOPE. The amplitude can be then better treated within the 
T matrix approach which sums up the ladder diagrams in the particle-particle chan- 
nel. Calculations of MNB processes with the vacuum T matrix [65] found essentially 
smaller emissivity then that given by FOPE. Also as has been demonstrated in [66] 
the in-medium modifications of the T matrix additionally suppress the rates of both 
MMU and MNB processes. Thus even at small densities the FOPE model may give 
only a rough estimate of the emissivity of two nucleon processes. At p > 0.5 — 0.7 p 
reactions in particle-hole channel and more specifically with participation of the soft 
pion mode begin to dominate. 

Evaluations [23,16] showed that the MMU is dominated by the second diagram of 
Fig. 1, namely by contributions to it given by the first two diagrams of the series of 
available diagrams 
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V 




(27) 



whereas the third diagram which naturally generalizes the corresponding MU(FOPE) 
contribution gives only a small correction for p > p 0l [23,42,43,16]. The resulting 
emissivity given by two first diagrams in a simplified notation reads as follows, see 

([16,26]), 



MMU 



[MOPE] ~ 2.4 x 10 24 T t 



10/3 (OH 



,Po. 



m 



N 



m-n 



acu 7r o(p Fn ) 



m-n 



«^tt± (PFn) 



x T 8 F X C(A„) C(A P ) 



erg 



cm 3 sec 



(28) 



where T 9 = T/10 9 K is the temperature, m* n and m* are the nonrelativistic effective 
neutron and proton masses, respectively, 



T 8 = 7Jj(cd ^PFe,q^PFeh 2 (9nn,U ^ 0, k = PFnh^id', 0, PFn) , 

^^ = TTM ' 

and the second term in the factor 



(29) 



2/3 



4>y 2 (g',0,p Fn ) \po 



(30) 



is the contribution of the pion decay from intermediate states (first diagram). The 
quantity T effectively accounts for a product of the NN correlation factors in dif- 
ferent ttNiN 2 vertices. For charged pions the value p n ^ is incorporated in the 
expression for the effective pion gap, for neutral pions p n = 0. The value of the 
parameter a is given by a = 1 for p < p C7r , and a = y/2 for p > p cn taking account 
of the new excitations on the ground of the tt condensate vacuum. The factor 



C(Aiv) 



exp {-A N /T) T < T cN , 
1 T>T cN , N 



{n,p), 



(31) 



estimates the suppression caused by the neutron-neutron and proton-proton pair- 
ings. A deviation of these factors from simple exponents can be incorporated as in 
ref. [33]. 
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The ratio of the emissivities of MMU(MOPE) and MU(FOPE) is roughly as follows 

"" 10 ^j-, r^j— t r \PIPo) • (32) 



rs_> 



e Mt/[ F0 PE] ^APFnH±(PFn) 

For p — po this ratio is ~ 10 whereas being estimated with the only third diagram 
(27) it would be less then unit. 

3.2 Emissivity of DU-like processes 



NPBF processes. The one-nucleon processes with neutral currents given by the 
second diagram (19) for N = n and p are forbidden at T > T c by energy-momentum 
conservation law but they can occur at T < T c . Then physically the processes relate 
to NPBF, see Fig. 2. However they need a special technique to be calculated, see 

[32,24]. 

These processes n — > nvv and p — > pvv play very important part in the cooling 
of superfluid neutron stars, see [24,25,16,26,34,33]. The emissivity for 3 types of 
neutrinos is given by [24,25] 3 



e[nPBF]n = 3 \ ^ (£i7 2 (/»n) + 6<&7 2 (g»n))Pf»<A£ 

~ C x 10 28 (p/p ) 1/3 (m;/m jV )(A„/Mey) 7 J(A„/r), f 8 , T < T cn , (33) 

cm d • sec 

where G ~ 1.17 -10~ 7 GeV~ 2 is Fermi constant of weak interaction, A„ is the neutron 
gap, T cn is critical temperature for neutron pairing, £i = 1, £2 = for 5-pairing and 
£1 = 2/3, £2 = 4/3 for P-pairing, cf. [33]. For typical neutrino energies uj q = q 
correlations are not so essential as it is for lu <$; q. Taking 7 2 ~ 0.9 in the range of 
5-pairing we get £ ~ 3 whereas with 7 2 ~ 0.7 for a more dense region of P-pairing 
we obtain £ ~ 6.6, in agreement with numerical evaluations [24] used within the 
cooling code in ref. [26], I{x) = J °° ch. 5 ydy/ (exp(rcchy) + 1) . One should notice 
that this numerical estimate is by a factor of (20 -r- 50) larger than that in [32]. The 
emissivity of the process p — >■ pvv is given by [24] 



e(p ■, w - v) = 3 X ^^ + ?0l + ^^ ;[A,/n T<T V , (34) 

Ap is the proton gap, T cp is critical temperature for pairing of protons, which are 
paired in 5-state in neutron star matter, The vertex corrections take into account 
nucleon-nucleon correlations, £3^1 is due to the second diagram (20) and has a 



Here I remove some misprints existed in [32,24,25,16]. 
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complicated structure [54,55]. In the process (34) a part of NN correlations is es- 
pecially important. One has K 2 p ~ 1 instead of a small c v ~ 0.006 factor in the ab- 
sence of the correlation contribution. Thereby, as has been originally demonstrated 
in [24,26], the emissivity of the process p — >■ pvv can be compatible with that for 
n — >■ nvv in dependence on the relation between A p and A n . 

Inclusion of superfluid effects for rather low temperatures does not change our gen- 
eral conclusion. But very efficient processes become to be NPBF processes which 
compete with MMU processes. The former win the content for not too massive stars 
in agreement with the prediction of [24]. 

Pion (kaon) condensate processes. The pion condensate can be of three types: 
71"+, 7!"+, and 7T° with different values of the critical densities p C7T = (p C7r ±, p C7r +, p C7r o, 
see [67]. Thus above the threshold density for the pion condensation of the given 
type, the neutrino emissivity of the MMU process (28) is to be supplemented by the 
corresponding PU processes 




(35) 



To be specific, considering charged pion condensation we use a simplified expression 
including the nucleon-nucleon correlation effect in the ttNN vertices [22,16], 



e[PU]*1.5 x 10 27 ^^^( PFe , PFe ) f(g'Ap Fn )T« sin 2 -^-. (36) 
m v m l N p cnr sec 

Here p > p cn and sin# ~ J2\uj 2 \/m1 for 6 <C 1. Of the same order of magnitude 
are emissivities of other possible tt condensate reactions. Since 7r +, ~ condensation 
probably reduces the energy gaps of the superfluid states by an order of magnitude, 
see [68], we may assume that superfluidity vanishes above p C7r . Finally we note that 
though the PU processes have genuinely one-nucleon phase-space volumes, their 
contribution to the resulting emissivity is suppressed relative to the DU by an ad- 
ditional 7 2 (#', 0,pFn) suppression factor due to existence of the extra (irNN) vertex 
in the former case. 

Fig. 4 compares the mass dependence of the neutrino cooling rates L v jCv, L u is 
neutrino luminocity, Cy is the heat capacity, associated with MMU (MOPE) and 
MU (FOPE) for non-superfluid matter. For the solid curves, the neutrino emissivity 
in pion-condensed matter (to be conservative we took threshold density to be 3p ) 
is taken into account according to eq. (36) and (28) with the parameter a = y/2. 
Dashed curves correspond to the case where no pion condensation is allowed. As one 
sees, the medium polarization effects included in MMU may result in three order 
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> FOPE with n-cond. 
i MOPE with 7r-cond. 
i MOPE without ji-cond. 

> FOPE without 7t-cond. 




Fig. 4. Cooling rate due to neutrino emission as a function of star mass for a representative 
temperature of T = 3 x 10 8 K. The solid curves refer to cooling via MU(FOPE) and 
MMU(MOPE). For p > 3/?o PU is included. The dashed curves refer to cooling without 
PU. Superfluidity is neglected. 

of magnitude increase of the cooling rate for the most massive stars. Even for stars 
of a rather low mass the cooling rate of MMU is still few times larger than for MU 
because even in this case the more efficient rate is given by the reactions shown 
by the right diagram in Fig.l (first two diagrams of eq. (27)). The cooling rates of 
the 1.8M mass models with and without pion condensate differ only by a factor 
of 5 in this model. If we would take p C7T to be smaller the ratio of emissivity PU to 
MMU would decrease and could even become ^ 1 in the vicinity of p C7T . One sees 
that the reaction rates for the FOPE model are rather independent of the star's 
mass for the stars with masses below the critical value 1.63Mq, at which transition 
into the pion condensed phase occurs, and show a flattening behaviour above the 
critical mass. It is to be stressed that incorporation of pion condensation processes 
seems quite inconsistent, whether one applies FOPE model since pion condensation 
arises at p > p C7T only as the consequence of the pion softening existed for p < p C7T . 
Opposite, the MOPE model [23] consistently takes into account the pion softening 
effects for p < p C7r and both the pion condensation and pion softening effects on the 
ground of the condensate for p > p C7r . For p > p cK the kaon condensate processes 
come into play. Although kaon condensation can be of both S and P wave types 
[69], most popular is the idea of the S wave K~ condensation (e.g. see [13]) which is 
allowed at the condition p, e > m* K - due to the possibility of the reaction e — >■ K~v. 
The corresponding condition for the pions \i e > m*_ is not fulfilled due to a strong 
S wave repulsion in the latter case. The neutrino emissivity of the K condensate 
processes (see first diagram of (35)) is given by equation analogous to (36) with a 
different NN correlation factor and an additional suppression factor due to a small 
contribution of the Kabibbo angle. However qualitatively scenarium that permits 
pion condensate processes is analogous for that with the kaon condensate processes. 

Other resonance processes. There are many other reaction channels allowed in 
the medium. E.g any Fermi liquid allows for propagation of zero sound excitations 
of different symmetry related to the pion and the quanta of a more local interaction 
determined via f a ^ and g a ^. These excitations being present at T ^ may par- 
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ticipate in the neutrino reactions. The most essential contribution comes from the 
neutral current processes [23] given by first two diagrams of the series 

v 
v v ^- 

V n » AL » ■ n 



o 



n l 



" ' + ... (37) 



Here the dotted line is zero sound quantum of appropriate symmetry. These are 
the resonance processes of DU-type analogous to those processes going on the con- 
densates with the only difference that the rates of reactions with zero sounds are 
proportional to the thermal occupations of the corresponding spectrum branches 
whereas the rates of the condensate processes are proportional to the modulous 
squared condensate mean field. Contribution of the resonance reactions is however 
rather small due to a small phase space volume (q ~ T) associated with zero sounds. 
Please also bear in mind an analogy with the processes (37) with the phonon pro- 
cesses in the crust. 

DU processes. The proper DU processes in matter, asn^ peTv e an d v e ~~ *• nv e 



K 




(38) 



P P 

should also be treated with the full vertices. They are forbidden up to the density p c u 
when triangle inequality pp n < Pf p + PFe begins to fulfill. For traditional equations 
of state like that given by the variational theory [18] DU processes are permitted 
for p > 5po- The emissivity of the DU processes renders 

e ov ^ 12 . 1027 <!p; f^—) 7 > in [ C (A n ), C(A P )] T 9 6 -^i-, (39) 

mjf VlOOMeV/ p cnr sec 

where pi = fi e = p^ is the chemical potential of the leptons in MeV. In difference 
with the usually exploited result [15], eq. (39) contains a pre-factor 7 | [26] due to 
taking into account of the NN correlations in the j3 decay vertices. 

It was realized in ref. [70] that the softening of the pion mode in dense neutron 
matter could give rise to a rearrangement of single-particle degrees of freedom with 
an abnormal increase of the effective non-relativistic masses in a density interval 
p C F < p < pc-K that in its turn, opens a DU channel of neutrino cooling of neutron 
stars from the corresponding layer. 

A rearrangement of single-particle degrees of freedom occurs if the necessary con- 
dition for stability of the normal state of a Fermi liquid is violated. At T = this 
condition states that the change of the ground-state energy must remain positive for 
any admissible variation Sn(p) of the Landau quasiparticle distribution n(p) away 
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from its normal-state form 9(p — pp), i.e., SE = / £ (p, n(p)) 8n{p)S^ > , where 
£(p, n(p)) = £ (p, n(p)) — /i is the energy of a quasiparticle, measured from the chem- 
ical potential p. This stability condition fails if a depression with £ < forms in 
the spectrum £(p) at p > pp ; it likewise fails if there arises an elevation with £ > 
at p < pp. The rearrangement is precipitated when the density p reaches a critical 
value p c p where there emerges a new root p = p of the equation £ (p, n(p); p c p) = , 
a relation that ordinarily serves merely specify the Fermi momentum pp. 

The nucleon quasiparticle spectrum is determined from the corresponding Dyson 
equation which in our MOPE model renders 



e(p) = ?(P) + A 2 / — — % ,^ 2/A , 2 n(k)-^ , (40) 

J u} 2 + -i m [{p-k) 2 -kl)Y/Akl (2tt)^ 



m 



where A 2 = /^v7 2 (<?mn 0; ko)m* N /m,]s[, and we expanded the value ui 2 (k) near ko, 7 
is the corresponding coefficient. Numerical study of [70] found the bifurcation point 
and documented three features of the rearrangement of the single-particle degrees 
of freedom, (i) The critical density p c p for the rearrangement is less than the critical 
density p C7T for pion condensation. Since both phenomena stem from the strong 
momentum dependence of the Landau amplitude, rearrangement of the quasiparticle 
distribution may be regarded as a precursor of pion condensation, (ii) The bifurcation 
point is positioned at p < 0.2p F for the density interval under consideration, (iii) 
The ratios p c /p C7r and Po/pf are insensitive to the actual value taken by p C7r within the 
current range of theoretical predictions, (iiii) For the momenta above the bifurcation 
point p > po there arises a dip in the profile £(p) with a minimum at p ~ 0.bpp n . 

An equation analogous to (40) can be applied to the proton subsystem of neutron- 
star matter. One is then dealing with the charged pion mode, and the neutrons form 
a strongly momentum-dependent external field in which the protons seek to occupy 
states with the lowest single-particle energies. The critical densities for neutral and 
charged pion condensation do not differ greatly, so we can appropriate the results 
obtained for neutron matter to argue that protons will leave the old Fermi sphere 
and occupy states of relatively large momentum in the region p ~ 0.bpp n . We defer 
here the impact of this further rearrangement. 

Beyond the bifurcation point, the creation of a bubble at low neutron momenta in the 
Fermi sea means that the triangle inequalities are satisfied without the conventional 
requirement that the proton fraction exceeds some 11-14% , cf. [15]. Accordingly, the 
triangle condition for momentum conservation is readily met and the DU processes 
become active in the density regime just short of the threshold for pion condensation. 
The neutrino emissivity goes, as usual, like T 6 . Therefore the rearrangement of the 
neutron Fermi sphere that we have envisioned serves to turn on the DU process at 
a lower density. 

There also exists a more radical scenario for rearrangement of the quasiparticle distri- 
bution [71-73]. In this case the condition £ = 0, that has determined the bifurcation 
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point continues to fulfill in a density interval p c p < p < p C7T and in the corresponding 
finite momentum interval pi < p < pf. In this interval the occupancy n(p) becomes 
partial, between and 1 that determines the effective nonrelativistic mass relation 
to the temperature as m* oc m^ep/T. The family of quasiparticles having momenta 
p e [Pi,Pf] characterizing by dispersionless spectrum (at T = 0) has been called 
the fermion condensate because of a conspicuous analogy with the low-temperature 
Bose gas, where the energy of condensate particles is also equal to the chemical po- 
tential fj,. It is worth noting that a key signature of fermion condensation has been 
observed in strongly correlated electron systems. Flat portions of the single-particle 
spectrum have been seen experimentally in a number of high-temperature supercon- 
ductors ([74]), and this phenomenon has been linked to fermion condensation ([72]). 
Numerical studies [75] demonstrate that fermion condensate wins over the burble 
formation at T ^ 0. Although the rearranged Fermi condensate momentum distri- 
bution differs from that of the bubble rearrangement, its structure will also permit 
thermal excitations to occur at low neutron momenta. Hence we may again expect 
most neutron stars to contain a region of relatively moderate density, bounded below 
by p C F and above by p C7r , in which the DU process operates vigorously. However, due 
to the new feature of a temperature-dependent neutron effective mass, M n oc 1/T, 
we may anticipate an extra enhancement of the neutrino emissivity relative to the 
standard result of [15], corresponding to a reduction in the power of the tempera- 
ture dependence from T 6 to T 5 . At early hot stage this layer becomes to be opacue 
to the neutrinos slowing the transport from the massive neutron star core to the 
exterior. 



3. 3 Comparison with soft X ray data 



At first seconds-hours neutrino opacity effects become to be not efficient and the 
neutron star is transparent for neutrinos/antineutrinos. The heat transport within 
the crust establishes homogeneous density profile at times & lOyr. After that time 
the subsequent cooling is determined by simple relation CyT = —L where Cy = 
J2i Cy,i and L = J2iLi are the sums of the corresponding partial contributions to 
the heat capacity (specific heat integrated over the volume) and the luminosity 
(emissivity integrated over the volume). 

The nucleon pairing gaps are rather purely known. Therefore one may vary them. 
The "standard" and "nonstandard" scenarios of the cooling of neutron stars of sev- 
eral selected masses for suppressed gaps are demonstrated in the left panel of Fig. 5, 
[26]. Depending on the star's mass, the resulting photon luminosities are basically 
either too high or too low to account for the bulk of observed pulsar luminosities. 
Situation changes, if the MMU-process with the medium modifications of the pion 
exchange interaction, appropriate vertex corrections, and the radiation from inter- 
mediate reaction states [23,16] is included. Now, the cooling rates vary smoothly 
with the star mass (see right panel of Fig. 5) such that the gap between standard 
and enhanced cooling is washed out. More quantitatively, by means of varying the 
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Fig. 5. Cooling of non-superfluid neutron star models of different masses constructed 
for the HV EOS. The two graphs refer to cooling via MU(FOPE) +PU (left) and 
MMU(MOPE)+PU (right). In both cases, pion condensation is taken into account for 
the solid curves where p > p cv chosen to be 3/?o- The dashed curves in the right graph re- 
fer to the LV 2 (pF n ) parametrization without pion condensation. The observed luminosities 
are labeled by dots. Possibility of Fermi sea rearrangement is suppressed. 

star's mass between 1 — 1.6 Mq, one achieves an agreement with a large number 
of observed data points. This is true for both choices of the ui 2 parameterization, 
independently whether pion (kaon) condensation can occur or not. Two parameter- 
izations differ only in the range which is covered by the cooling curves. The only 
point which does not agree with the cooling curves belongs to the hottest pulsar 
PSR 1951+32. Other three points which according to Fig. 5, right panel, are also 
not fitted by the curves can be, nevertheless, fitted introducing slight changes of 
the parameters of the model. The high luminosity of PSR 1951+32 may be due to 
internal heating processes, which lead to a delayed cooling for star ages t > 10 5 yr, 
cf. [10]. 

We turn now to cooling simulations where the MU, NPBF, DU and PU take place 
simultaneously. Parameters of the pairing gaps are from Fig. 6 of ref. [26]. Fig. 6 
shows the cooling tracks of stars of different masses, computed for the HV equation of 
state. Possibility of Fermi sea rearrangement is suppressed. Very efficient at T < T c 
become to be NPBF processes which compete with MMU processes. The former win 
the content for not too massive stars in agreement with prediction of [24]. The DU 
process is taken into account in the right graph, whereas it is neglected in the left 
graph. The solid curves refer again to the o> 2 -parameterization with phase transition 
to a pion condensate, the dashed curves to the one without phase transition. For 
masses in the range between 1.0 and 1.6 M , the cooling curves pass through most 
of the data points. We again recognize a photon luminosity drop by more than two 
orders of magnitude of the 1.7 M mass star with pion condensate, due to vanishing 
of the superfluidity as a consequence of pion condensation. This drop is even larger 
if the DU is taken into account (right graph, this allows to account for the photon 
luminosity of PSR 1929+10). 
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Fig. 6. Cooling of neutron stars with different masses constructed for the HV EOS. The 
cooling processes are MU-VS86, PU (only solid curves), NPBF, PPBF, and DU (only in 
the right graph). The dashed curves refer to the M = 1.7 and 1.9M models without pion 
condensate. 

A comparison with the observed luminosities shows that one gets quite good agree- 
ment between theory and observation if one includes into consideration all awailable 
in-medium effects assuming that the masses of some of the underlying pulsars are 
different from the canonical value, M = 1.4 M Q . Since the masses of these pulsars 
are not known, no further conclusions about the actually operating cooling mecha- 
nism of these pulsars can be drawn yet. Future mass determinations of these objects 
will change the situation. 

Thus, ref. [26] demonstrated that the data can be treated with the help of the model 
based on the results which incorporate most important in-medium effects. We also 
point out that the description of these effects is constructed in essentially the same 
manner for all the hadron systems as neutron stars, atomic nuclei and heavy ion 
collisions, cf. [16]. 



3.4 Neutrino opacity 



Important part of in-medium effects for description of neutrino transport at initial 
stage of neutron star cooling was discussed in [23,24,16], where correlation effects, 
pion softening and pion condensation, the latter for p > p c , were taken into account. 
The neutrino and antineutrino mean free paths can be evaluated from the corre- 
sponding kinetic equations via the neutrino/antineutrino widths ^ v (&) = — 2ImII^ P N, 
where 11^ is the retarded selfenergy, or within the quasiparticle approximation for 
the nucleons they can be calculated via the squared matrix elements of the cor- 
responding reactions. Thereby the processes which the most efficiently contribute 
to the emissivity (at T ~ lMeV typical for the value of T opac ) are at early times 
responsible for the opacity. 
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In above scenarium at T > T c the most essential contribution was given by MMU. 
Taking into account of nucleon-nucleon correlations in the vertices of two-nucleon 
processes like MMU and MNB suppresses the rates, whereas account of the pion 
softening in the pion propagator leads to their essential enhancement. For rather 
massive neutron stars MOPE wins the competition. The mean free path of neu- 
trino/antineurtino in MMU processes is determined by the same diagrams (27) as 
the emissivity. Its calculation with the two first diagrams gives the following evalu- 
ation (see (28)) 



x mmu 15 . 10 5 / Po y 10 / 3 „4 
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From the relation X v ~ R follows evaluation of T opac . With the only first diagram 
we get a simple relation 



-opac _ ^P0^ 2 (p 2 Fn) m N 



15- : 7^^^f (42) 
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For averaged value of the density p ~ p corresponding to medium heavy neutron 
stars (< 1.4M ) with Cj 2 {p Fn ) ~ 0.8ml, r - °- 4 we S et T o P ac ^ 1.5MeV that is still 
smaller then the corresponding value 22MeV estimated with FOPE. For p ~ 2p 
that corresponds to a more massive neutron star we evaluate T opac ~ 0.5 MeV. Thus 
pion softening results in a substantial decrease of neutrino path lengths for typical 
energy ~ T and the value of T opac . 

The diffusion equation determins the characteristic time scale for the heat transport 
of neutrinos from the hot zone to the star surface t ~ R 2 Cv<7~ 1 T~ 3 /X 1/ , whereas 
it follows that t -10 min. for T ~ lOMeV and p ~ p , u 2 ~ 0.8m,,, T ~ 0.4, 
m* N /rriN ~ 0.9 and t — becomes as large as several hours for p ~ 2 -i- 3p . These 
estimates demonstrate that more massive neutron stars cool down more slowely at 
T > T opac and faster at subsequent times then the less massive stars. 

Due to in-medium effects neutrino scattering cross sections on the neutrons shown 
by the first diagram (19) require the factor 



2/ x _ l 2 (fnn,u,q) + 3g 2 A 7 2 (g nn ,u),q) 
7 W K Q) ~ Y^g\ ' (43) 

due to nucleon-nucleon correlations, as follows from eq. (17). This results in the 
corresponding suppression factor 0.5 -r- 0.7 of the cross sections. Neutrino scattering 
cross sections on the protons are modified differently, by 



K 2 (f np , fnn, uj , q) + 3g 2 A j 2 p {g nn , u, q) 



llM Q) = ~^^"^-\ W ^ (44) 
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that leads to analogous estimate of their contribution to the cross sections of the 
up — >■ up processes. 

Also there is an additional suppression of the uN scattering and MNB reaction rates 
for soft neutrinos {uj ^ 3T) due to multiple NN collisions as 



I? 1 



^Z (45) 



(Landau-Pomeranchuk-Migdal effect). One may qualitatively estimate these ef- 
fects multiplying the corresponding squared matrix elements by suppression factor 
Co(lo) = uj 2 /[uj 2 -\-T 2 n ] where T N is the nucleon width and the uj is the energy of uu or 
uu pair. However in general the prefactors are different for different diagrams under 
consideration. They follow from derivation from the general result [28]. Factor like 
C was used in refs [78,79] being introduced there at ansatz level. One does not need 
any ansatze reductions. Optical theorem formalism suggested in, see [28], allows to 
calculate the rates using the exact sum rule! The modification of the charged current 
processes is unimportant since the corresponding uj is ~ pp e 3> Fa- . 

The main physical result is that besides a suppression due to nucleon correlations 
and infra-red pre-factors, there exists an enhancement due to pion softening. The 
latter demonstrates that for densities p > p c \ ~ (0.5 — 0.7)p nucleon system begins 
to feel that it may have tt condensate phase transition for p > p c , although this p c 
value might be essentially larger than the value p c \. 



4 How to calculate the rate of radiation from dense medium 



Quasiparticle approximation for fermions is applicable if the fermion width is much 
less than all the typical energy scales essential in problem under consideration 
T^ <$; u c h- In calculation of the emissivities of uu reactions the minimal scale is 
uj ch ~ 3T, averaged energy typical for MNB reactions. For MMU it is ~ ppe. For 
radiation of soft quanta of fixed energy uj < T it is uj. Within quasiparticle approx- 
imation for fermions the reaction rate with participation of in-medium fermion and 
in-medium boson is given by [22,23] 



< 




> (46) 



For equilibrium (T 7^ 0) system there is the strict relation 

(< i&ipi >)(p) = iD- + = J^—, T B = -2ImD R , (47) 
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where (< (p\(f>\ >){p) means the Fourier component of the corresponding non- 
equilibrium Green function. Thus the rate of the reaction is related to the boson 
width being determined by the corresponding Dyson equation for D~ + having the 
form (9). The boson width T B is non-zero at the spectrum branches related to res- 
onance processes, like zero sound. The poles related to upper branches do not con- 
tribute at small temperatures due to a tiny thermal population of those branches. 
There is also a contribution to ImD^ proportional to Imll^ given by the particle- 
hole diagram. Taking Im part means the cut of the diagram. Thus we show [23] 
that this contribution is the same as that could be calculated with the help of the 
squared matrix element of the two-nucleon process 

<*£- 2 

(48) 



This is precisely what one could expect using the optical theorem. Thus unlimited 
series of all possible diagrams with in-medium Green functions (see (19), (35), (37), 
(38)) together with two-fermion diagrams (as given by (27)) and multiple- fermion 
diagrams (like (45)) would lead us to a double counting. The reason is that per- 
mitting the boson width effects (and beyond the quasiparticle approximation for 
fermions also permitting finite fermion width) the difference between one-fermion, 
two-fermion and multiple-fermion processes in medium is smeared out. All the states 
can still participate in subsequent production and absorption processes. Within 
quasiparticle approximation for fermions these further processes are suppressed for 
fermions (due to assumption IV "C e c /i) but allowed for the bosons. Staying with 
the quasiparticle picture for fermions, the easiest way to avoid mentioned double 
counting is just to calculate the reaction rates according to (46) thus with the help 
of the diagrams of the DU-like type, which already include all the contributions of 
the two-nucleon origin. On other hand, it is rather inconvenient to explicitly treat 
all one-nucleon processes dealing with different specific quanta instead of a more 
convenient treatment of the full NN interaction like in (48) which has a compli- 
cated structure. Besides, as we have mentioned, consideration of open fermion legs 
is only possible within the quasiparticle approximation for fermions since Feynman 
technique is not applicable anymore if Green functions corresponding to ingoing and 
outgoing fermions have widths that in other language means possibility of additional 
processes. Thus the idea came [24,28] to integrate over all in-medium states allowing 
all possible processes instead of specifying different special reaction channels. 

In refs [24,28,27] it was shown that OTF in terms of non-equilibrium Green's func- 
tions is efficient tool to calculate the reaction rates including finite particle widths 
and other in-medium effects. Applying this approach to the anti-neutrino-lepton 
(electron, \i~ meson, or neutrino) production [24] we can express the transition 
probability in a direct reaction in terms of the evolution operator S 

^ = *Js^e<«i*i« + *><« + .vi*io>. m 
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where we write explicitly the phase-space volume of vl states; lepton occupations 
of given spin, ni, should be put zero for v which are supposed to be radiated di- 
rectly from the system. The bar denotes statistical averaging. The summation goes 
over complete set of all possible intermediate states {X} constrained by the energy- 
momentum conservation law. It was also supposed that electrons (and muons) can 
be treated in the quasiparticle approximation, i.e. with zero widths. Thus, although 
possible, there is no need to consider them in intermediate reaction states. Mak- 
ing use of the smallness of the Fermi weak-interaction constant G, we expand the 
evolution operator as S ~ 1 — i f*™ T {Vw(x) S nuc i(x)}dx , where Vw is the Hamil- 
tonian of the weak interaction, in the interaction representation, S QUC \ is the part 
of the S matrix corresponding to the nuclear interaction, and T{...} stands for the 
chronological ordering operator. After substitution of this S matrix into Eq. (49) and 
averaging over the arbitrary non-equilibrium state of a nuclear system, there appear 
chronologically ordered (G ), anti-chronologically ordered (G ++ ) and disordered 
(G + ~ and G~ + ) exact Green functions. 

In graphical form the general expression for the probability of the lepton (electron, 
muon, neutrino) and anti-neutrino production is determined by the diagram 

I 




which represents the sum of all closed diagrams (—iU~ + ) containing at least one 
( — h) line. We will deal with full Green's functions and we will explicitly present 
all the contributions containing full G + ~ and G~ + Green's functions, cf. [28]. Then 
the contributions of specific processes contained in a closed diagram can be made 
visible by cutting the diagrams over the (H — ), ( — h) lines. In the framework of 
the quasiparticle approximation for the fermions (nucleons in our case) G~ + = 
—n5(e + \x — 6p — ReS B (e + /J,,p)) and the cut eliminating the energy integral thus 
has clear meaning. This way one establishes the correspondence between closed 
diagrams and usual Feynman amplitudes. In general case of finite fermion width the 
cut has only a symbolic meaning. 

Various contributions from {X} can be classified according to the number N of G~ + 
lines in the diagram. We can write, e.g., 



dWl? d 3 q p dqf 




(50) 



dt (2tv) 6 4:Uj j} uji 

v I v I 

This procedure suggested in [24] is actually very helpfull especially if the quasipar- 
ticle or quasiclassical approximations hold for the fermions (for loi + u 9 ~ T ^$> Fp , 
where T F ~ T 2 /e F is the fermion (nucleon) width). E.g. in the quasiparticle approx- 
imation any extra G~ + brings a small (T/e F ) 2 factor to the emissivity of the process. 
Dealing with small temperatures one can restrict by the diagrams of the lowest order 
in [G h G + ~) not forbidden by energy-momentum conservations putting T = in 
all G ++ and G Green functions. 
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Considering one-nucleon processes we decompose the first term in Eq. (??) as 
I v I v I v 

*ex- = xx + xx • <»> 

v I v I v I 

The full vertex in the diagram (51) of given sign is irreducible with respect to the 
(H — ) and ( — h) nucleon-nucleon-hole lines. This means it contains only the lines 

of a given sign, all ( ) or (++). Second diagram with anomalous Green functions 

exist only for systems with pairing. In the framework of quasiparticle approximation 
namely these diagrams determine the proper DU process and also NPBF processes 
calculated in [24,25] within this approach. In the quasiparticle picture the contri- 
bution to DU process corresponding to the one-loop diagram vanishes for p < p cU . 
Then the second term in (50) comes into play which within the same quasiparti- 
cle picture contains two-nucleon processes with one (G~ + G + ~) loop in intermediate 
states. 

The set of diagrams for n~ + can be further explicitly presented as series [28] 

+ -00 + -<x>- + Offi- * 



(52) 



where brief notations are used. Full dot is the vertex including all the diagrams of 
given sign, NN interaction block is a full block including diagrams of one sign. The 
lines are full Green functions with the widths. In general case given by (52), each 
diagram represents a whole class of perturbative diagrams of any order in the inter- 
action strength and in the number of loops and contains multiple-collision processes. 
The most essential term is the one-loop diagram, which is positive definite, and in- 
cluding the fermion width corresponds to the first term of the classical Langevin 
result, see [28]. The other diagrams represent interference terms due to rescatter- 
ings. In some simplified representations (being often used) the 4-point functions 
(blocks of NN interaction of given sign diagrams) behave like intermediate bosons 
(e.g. phonons or dressed pions). In general it is not necessary to consider different 
quanta dealing instead with the full A^A" interaction (all diagrams of given sign). 
For particle propagation in an external field, e.g. infinitely heavy scattering centers 
(proper Landau-Pomeranchuk-Migdal effect), only the one-loop diagram remains, 
since one deals then with a genuine one-body problem. In the quasiparticle limit 
the diagrams 1, 2, 4 and 5 of (52) correspond to the MMU and MNB processes 
related to a single in-medium scattering of two fermionic quasi-particles and can be 
symbolically expressed as Feynman amplitude (53a) 
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(c) 
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The full circle denotes the effective weak coupling vertex. The one-loop diagram in 
(52) is particular, since its quasiparticle approximant in many cases vanishes as we 
have mentioned. However the full one-loop includes quasiparticle graphs of the type 
(53b), which survive to the same order in F/e as the other diagrams (therefore in 
[24] where quasiparticle picture was used this diagram was considered as allowed 
diagram). In quasiparticle series such a term should be included in second diagram 
of (refqp) although beyond the quasiparticle approximation it should be included 
as self-energy insertion to the one-loop result [28]. In fact it is positive definite and 
corresponds to the absolute square of the amplitude (53a)). The other diagrams 2, 
4 and 5 of (52) describe the interference of amplitude (53a) either with those where 
the weak coupling quantum couples to another leg or with one of the exchange 
diagrams. Thereby for neutral interactions diagram (52:2) is more important than 
diagram 4 , while this behavior reverses for charge exchange interactions (the latter 
is also important for gluon radiation from quarks in QCD transport due to color 
exchange interactions) . Diagrams like 3 describe the interference terms due to further 
rescatterings of the source fermion with others as shown by (53c). Diagram (52:6) 
describes the production from intermediate states and is given by Feynman graph 
(53d). For photons in the soft limit (iO q <C e F ) this diagram (53d) gives a smaller 
contribution to the photon production rate than the diagram (53a) in quasiparticle 
approximation, where the normal bremsstrahlung contribution diverges like l/u) q 
compared to the 1/e^-value typical for the coupling to intermediate fermion lines 
[23]. For vv bremstrahlung it gives zero. However in some specific cases the process 
(53d) might be very important even in the soft limit. This is indeed the case for 
MMU process which is of prime importance in the problem of neutrino radiation 
from the dense neutron star interior. Some of the diagrams which are not presented 
explicitly in eq.(52) give more than two pieces, if being cut, so they do not reduce 
to the Feynman amplitudes. 

With Fp ~ n 2 T 2 /e F for Fermi liquids, the quasiparticle approximation for fermions 
constitutes a consistent scheme for all thermal excitations Ae ~ T <$; e F . However 
with the application to radiation at somewhat higher temperatures this concept is 
no longer justified. In particular, the remaining series of quasiparticle-diagrams is 
no longer convergent unless u > F F , since arbitrary powers in F F /u appear, and 
there is no hope to ever recover a reliable result by a finite number of quasiparticle 
diagrams for the production of soft quanta! With full Green's functions, however, 
one obtains a description that uniformly covers both the soft (u <^i F F ) and the 
hard (lu 3> IV) regime. 

In order to correct quasiparticle evaluations of different diagrams by the fermion 
width effects for soft radiating quanta one can simply multiply the quasiparticle 
results by different pre-factors [28]. E.g., comparing the one-loop result at non-zero 
T^ with the first non-zero diagram in the quasiparticle approximation (IV = in 
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the fermion Green's functions) we get 




(54) 



QPA 



at small momentum k where Co is given by (??) which cures the defect of the QPA 
for soft co. This factor complies with the replacement u — > co + ir. A similar factor 
has been observed in the diffusion result, where however the macroscopic relaxation 
rate T x enters, due to the resummation of all rescattering processes. 

Along similar routes the correction factors for the higher order diagrams can be 
derived. For the next order diagrams we have 
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In the quasiclassical limit for fermions (with occupations rip ~ exp((e — np)/T)) all 
the diagrams given by first line of series (52) with an arbitrary number of — h NN- 
interaction insertions can be summed up leading to the diffusion result. For small 
momenta q this leads to a suppression factor of the form C = co 2 /(co 2 + T 2 ,). In 
general case the total radiation rate result is obtained by summation of all diagrams 
in (52). 

Concluding, we demonstrated how to calculate the rates of different reactions in 
dense equilibrium (and non-equilibrium) matter and compared the results derived 
in closed diagram technique with those obtained in the standard technique of com- 
puting of the squared matrix elements. The nuclear medium cooling scenario allows 
easily to achieve agreement with existing data. However there remains essential 
uncertainty in quantitative predictions due to a pure knowledge of the residual in- 
teraction above treated within a phenomenological Fermi liquid model. 
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